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A.  SCOPE 

BI0DYN-80  is  a  versatile  computational  tool  used  to  determine  trans- 
missibilities  (transfer  functions)  between  vertical  and/or  fore-aft 
vibrational  inputs  and  important  biodynamic  outputs,  such  as  motions  of 
the  torso,  head,  eyes,  arms  or  hands.  The  program  scenario  assumes  a 
seated  pilot,  gripping  an  arbitrary- angle  stick  and  viewing  a  display, 
possibly  engaged  in  a  tracking  task.  The  physical  model  uses  an  "iso¬ 
morphic,"  (lumped  parameter,)  approach  to  represent  the  relevant  por¬ 
tions  of  the  whole-body  torso,  limbs  and  head,  as  well  as  postural 
compliances  among  the  joints.  The  implementation  of  this  model  includes 
a  chain  of  Interacting  parallel  and  serial  second-order  elements,  with 
neuromuscular  and  other  force  feedbacks  at  the  arm  or  head.  The  result¬ 
ing  equations  are  in  "second-order  element"  matrix  form  and  apply  to  a 
wide  range  of  seated  postures.  A  separate  input  file,  which  describes 
the  particular  set  of  parameters  to  be  used,  is  created  by  the  user 
(usually  by  modifying  one  of  a  cataloged  set).  This  file  is  incorpor¬ 
ated  in  the  matrix  to  produce  linearized  coefficients  for  perturbations 
about  the  selected  equilibrium  posture.  A  variety  of  outputs  and  inputs 
can  be  specified  to  evaluate  the  desired  transmissibility  transfer  func¬ 
tions,  and  these  are  written  to  a  file  in  formatted  form  for  plotting  or 
use  in  other  programs.  On-line,  printer-drawn  frequency  response  plots 
(Bode  format)  are  available  to  aid  in  interactive  user-computer  opera¬ 
tions,  or  to  screen  significant  results  from  batch  runs. 

■  .  BACXGEOOED 

BI0DYN-80  is  one  result  of  a  several  year  small-scale  development 
effort,  reported  in  detail  in  References  1-6.  It  is  based  on  vibration 


Measurements  made  at  the  Aeroapace  Medical  Research  Lab/Blodyaamics 
Division  (AMRL/BB)  and  elsewhere.  Most  of  the  torso,  limb  and  stick 
model  elements  are  based  on  independent  vibration  measurements  (e.g., 
Reference  5),  and  the  neck,  head  and  eye  effects  show  promising  corre¬ 
lations  with  the  few  available  measurements  on  Image  motion  effects 
(Reference  6).  However,  many  aspects  remain  to  be  explored,  validated 
or  upgraded  as  more  experiments  are  run  and  interpreted  via  BI0DYN-80. 

C.  APPLICATIONS 

The  possible  applications  of  this  program  are  many.  It  should  be 
used  in  the  early  stages  of  experimental  design  for  determining  the 
optimal  locations  for  vibration  measurements  and/or  selection  of  fre¬ 
quencies.  It  can  also  be  used  by  development  engineers  for  solving 
practical  pilot-vehicle  interface  design  problems  such  as  pilot-induced 
oscillations  and  for  optimizing  design  alternatives  such  as  seat  loca¬ 
tion,  orientation  and  suspension  parameters.  Flight  control  system 
designers  can  make  use  of  BIODYN-80  output  in  optimizing  vehicle/  air¬ 
crew  ride  qualities  and  visual  performance  effects,  possibly  incorporat¬ 
ing  anti-vibration  devices  to  improve  the  design. 

One  further  application  of  BIODYN-80  deserves  special  mention.  Its 
vibration-input  to  blodynamlc  parameter-output  transfer  functions  are 
ideally  suited  as  input  to  PIVIB,  another  software  package  which  relates 
pilot  tracking  performance  to  the  vibration  environment  (Reference  7). 
PIVIB  accepts  biomechanical  transfer  functions  in  the  format  created  by 
BIODYN-80.  The  details  of  the  BI0DYN-80/PIVIB  interface  will  be  found 
in  a  later  section. 


The  remainder  of  this  report  details  the  use  of  BIODYN-80,  and  des¬ 
cribes  the  model  and  the  equations  comprising  it.  Detailed  instructions 
for  the  creation  of  the  required  input  files,  and  a  complete  example 
problem  are  also  included. 
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SOFTWARE  OVERVIEW 


Figure  1  19  a  functional  block  diagram  description  of  the  ele¬ 

ments  In  BIODYN  and  Its  Interface  with  PIVIB.  The  BIODYN-80  package  Is 
composed  of  three  programs.  The  first,  called  CREATE,*  Interactively 
sets  up  the  two  input  files  used  by  BIODYN.  The  second  program  called 
BIODYN,  is  the  actual  "number  cruncher,"  which  structures  and  solves  the 
biomechanical  equations  and  computes  the  desired  transfer  functions. 
The  third  program  called  PLOT,  accepts  the  file  of  BIODYN  transfer  func¬ 
tions,  prints  selected  ones  in  a  form  readily  comprehended  by  the  user, 
and  prepares  "quick  plot"  Bode  plots  on  the  line  printer,  to  facilitate 
a  visual  interpretation  of  the  transfer  function  information.  Both 
BIODYN  and  PLOT  are  designed  primarily  as  batch  programs  while,  CREATE 
permits  conversational  user  interaction  in  structuring  input  data. 

A  subsequent  link  in  this  series  of  programs  is  PIVIB.  It  Is  a 
large  batch  program  with  three  modules.  The  first,  BDMOD ,  computes  the 
response  behavior  of  the  various  biomechanical  subsystems.  The  second, 
PVMOD ,  uses  the  results  of  BDMOD  and  the  BBN  optimal  control  model  to 
estimate  pilot  tracking  performance  within  the  vibration  environment. 
The  final  module,  VEXEC,  provides  the  top  level  communication  interface 
between  BDMOD  and  PVMOD,  and  performs  no  actual  computation.  The  BDMOD 
module  expects  biodynamic  transfer  functions  in  a  format  generated  by 
BIODYN.  Only  the  relationships  of  BIODYN-80  to  PIVIB  are  described 
herein,  because  PIVIB  is  run  separately  and  has  its  own  User's  Manual 
(Reference  7). 

Figure  2  presents  a  flow  chart  covering  the  use  of  BIODYN  and  PIVIB 
in  a  given  session.  The  detailed  description  of  the  various  steps  for 
BIODYN  is  found  in  Section  D.  Note  that  the  flow  of  execution  can  be 


*The  CREATE  program  should  not  be  confused  with  INTERCOM'S  EDITOR 
subcommand  CREATE;  it  has  been  given  the  permanent  file  name  EXECRT  to 
distinguish  it  from  the  latter. 
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used  to  solve  a  single  problem,  by  submitting  a  single  batch  request,  or 
to  Iterate  on  a  design  by  submitting  a  number  of  batch  requests  using 
the  same  CHOICES  file  and  slightly  altered  PARAMETER  files  (heavy  lines 
In  Figure  2). 

BIODYN-80  uses  three  separate  files  In  performing  its  computations. 
The  PARAMETER  file  contains  the  set  of  96  parameters  used  to  define  the 
specific  pllot/posture/dlsplay /vibration  characteristics.  Appendix  A 
contains  a  complete  description  of  each  parameter,  Including  definition, 
mnemonic,  nominal  value,  recommended  range  of  values  and  reference 
(where  available).  CREATE  Is  used  to  assemble  this  file  and  can  modify 
an  existing  file  or  produce  an  entirely  new  file.  The  CHOICES  file 
contains  the  list  of  desired  transfer  functions  to  be  computed  and 
output  by  BIODYN,  as  well  as  directives  for  producing  the  line  printer 
Bode  plots.  Again,  this  CHOICES  file  Is  assembled  by  CREATE.  Finally, 
the  TF  file  Is  used  to  store  the  resulting  transfer  functions  output  by 
BIODYN,  and  is  read  by  the  PLOT  routine  for  generating  Bode  plots. 

PIVIB  employs  a  single  large  file  to  direct  Its  flow  of  execution. 
This  file  defines  the  vibration  environment,  biomechanical  transfer 
functions,  tracking  dynamics,  tracking  performance  requirements,  and 
pilot  limitations  (bandwidths,  time  delays,  etc.).  Currently,  this  file 
Is  assembled  in  the  editor,  using  output  from  BIODYN-80  if  desired. 

BIODYN-80  and  each  of  its  predecessors  were  developed  on  the  Tym- 
share,  Inc.,  PDP-10  computers.  BIODYN-80  has  been  adapted  to  the  CDC 
"Intercom"  System  on  the  CDC  6600  or  CYBER  175  at  WPAFB  In  order  to 
increase  Its  availability  to  Air  Force  users,  and  to  Interact  with 
PIVIB,  which  is  also  operable  on  the  WPAFB  CDC  computer.  The  details  of 
this  manual  will  address  Its  use  on  the  CDC  machine;  an  example  of  a 
Tymshare  session  is  given  in  Appendix  D.  Throughout  the  manual,  how¬ 
ever,  it  Is  assumed  that  the  user  is  familiar  with  the  WPAFB  CYBER  175, 
and  in  particular  has  experience  with  the  INTERCOM  operating  system.  If 
not,  the  user  should  read  References  8  and  9  first. 
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MODEL  DESCRIPTION 


Three  distinct  subsystem  models  are  included  in  BI0DYN-80.  They  ere 
described  individually  below. 

A.  BIOMECHANICAL  MODEL 

Figure  3  (updated  from  Figure  2  of  Reference  3)  presents  the 
biomechanical  model  and  defines  many  of  the  necessary  parameters  that 
describe  the  nominal  (or  trim)  situation.  It  utilizes  an  "isomorphic, " 
or  body-mimicking  representation,  of  the  major  body  segments  in  their 
orientations,  simplified  to  a  minimum  number  of  lumped  parameter 
equivalents.  The  biomechanical  features  Include: 

•  Semisupine  torso;  sliding  hip  pluB  rocking  chest 
supported  on  a  compliant  buttocks/seat. 

•  Head  bobbing  on  an  articulated  neck  with  passive 
compliance,  or  active  neuromuscular  system. 

•  Upper  arm  and  forearm  links  plus  grip-interface  com¬ 
pliance,  driven  by  an  active  neuromuscular  system. 

•  Arm-rest  restraints  (optional). 

•  Stick  "feel  system"  dynamics  from  zero  to  infinite 
stiffness,  and  any  angle  of  stick  or  grip. 

The  simplified  torso  model  was  derived  to  describe  the  dominant  mo¬ 
tions  of  the  head  and  arm  elements;  the  "pin  joint"  node  between  upper 
and  lower  torso  segments  is  not  meant  to  represent  any  physical  feature. 
In  practice,  the  masses  and  inertias  are  obtained  from  tabulated  bio¬ 
mechanical  and  anthropometric  data  for  the  appropriate  sized  person 
(e.g.,  Reference  11),  the  postural  angles  are  based  on  the  actual 
situation  (preferably  via  a  side-view  photo),  and  the  spring  forces  and 
damping  coefficients  are  fitted  to  data  or  taken  from  other  sources 
(e.g..  References  3  and  12). 
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Figure  3  Main  Biomechanical  F.lements 


B.  LIMB  NEUROMUSCULAR  MODEL 


The  active  neuromuscular  system  noted  on  Figure  3  Is  a  schematic 
representation  of  the  net  effect  of  complex  agonlst/antagonlat  muscle 
pairs  controlling  the  upper  an  or  head  based  on  the  work  summarized  In 
Reference  10.  An  ”NM  switch"  Is  defined  In  the  PARAMETER  file  which 
causes  this  neuromuscular  model  to  control  the  limb  (NM  -  0)  or  the  head 
(NM  "l).  A  linearized  representation  of  the  limb  neuromuscular  model 
Is  shown  In  Figure  4,  while  the  head  neuromuscular  model  is  depicted  In 
Figure  5.  This  model  relates  the  action  of  the  muscle  pairs  to  the 
effective  (spindle)  sensors  of  muscle  length  and  force  as  well  as 
proprioception  from  the  stick  grip  interface  (in  the  case  of  the  limb 
neuromuscular  model)  or  the  head-neck  interface  (in  the  case  of  the  head 
neuromuscular  model),  thus  closing  the  receptor-CNS-effector  loop. 

Unless  neuromuscular  properties  are  being  Investigated,  It  Is  recom¬ 
mended  that  the  typical  values  of  the  parameters  shown  In  Table  A-l  be 
used.  These  are  representative  of  a  normal  person's  arm-hand  or  head- 
neck  system,  and  generally  yield  reasonably  damped  neuromuscular  servo 
properties.  The  neuromuscular  parameters  listed  In  Appendix  A  are 
characteristic  of  the  largest  muscles  in  the  body  (e.g.,  the  legs),  but 
experience  has  shown  that  the  dynamic  properties  (torque/lnertla  ratios, 
damping  ratios,  natural  frequencies,  etc.)  are  about  the  same  for  all 
postural  muscle  pairs.  Here,  an  empirical  scale  factor  S  (>  1.0)  is 
used  to  scale  the  normalized  muscle  to  a  particular  configuration,  as 
though  the  muscle  acted  normal  to  the  upper  arm  c.g.  A  more  detailed 
description  of  the  neuromuscular  model  is  found  in  Reference  10. 

C.  HECK/HEAD/EYE  MODELS 

Figure  6  illustrates  the  basic  elements  involved  in  the  neck/head/ 
eye  model.  Details  involving  validation  and  example  use  of  this  subsys¬ 
tem  model  are  found  in  Reference  6.  The  task  is  to  keep  the  eye  fixated 
on  the  target,  i.e.,  null  the  relative  eye  (point  of  regard)  deviation 
(RED)  at  the  display.  The  moving  base  can  produce  image  error  distur¬ 
bances,  both  from  induced  head  rotation  and  translation  as  well  as 
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Simplified  Head  Pointing  Neuromuscular  System 
(When  Option  NM  »  1) 


s)(2  -  rus)  |  Angular  Image  Motion  Feedback 


I 

i 
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Figure  6  Block  Diagram  for  Head,  F.ve,  Vestibular  Models  and  Display  Kinematics 


target  motions.  The  "vestibulo-ocular  reflex”  is  a  "crossfeed"  (or 
"feedforward”)  that  rotates  the  eye  oppositely  to  the  head  to  compensate 
for  head  rotation,  i.e.,  to  maintain  approximate  inertial  eye  fixation. 

The  "Fixation  Reflex"  or  feedback  tracking  loop  involves  the  sub¬ 
ject's  efforts  to  null  the  image  error  0j  by  compensatory  eye  movements. 
The  "Target  Pursuit"  path  (shown  dashed)  models  the  subject's  "feedfor¬ 
ward"  operations  on  the  perceived  absolute  target  motions  in  inertial 
space  (as  distinct  from  the  image  error  motions).  For  highly  predict¬ 
able  and  perceivable  target  motions  the  Target  Pursuit  path  is  capable 
of  greatly  reducing  the  image  errors,  but  it  is  not  involved  in  most 
vibration  cases  because  of  insufficient  bandwidth. 

The  linearized  model  formulation  assumes  that  the  target  remains 
within  the  foveal  area  (3-4  deg  field),  and  that  angular  velocities  are 
small  (i.e.,  less  than  20  deg/sec.)  such  that  saccades  are  rare  and  can 
be  ignored.  Additionally,  the  vibration-induced  motions  will  always  be 
small  enough  to  permit  linearization  of  all  angular  functions  and  allow 
a  quasi-linear  representation  of  the  dynamic  elements  about  each  "oper¬ 
ating  point"  (posture,  view  geometry,  frequency).  The  operating  point 
angles  can  have  any  values  within  human  limits  for  a  seated  pilot  or 
crewman. 

The  three  subsystem  models  described  above  contain  parameters  which 
must  be  assigned  values  via  the  PARAMETER  file,  in  order  to  describe  the 
desired  problem  to  BIODYN.  These  parameters  are  defined  in  Appendix  A. 
The  available  input  (forcing  function)  variables  and  output  (response) 
variables  are  given  in  Appendix  B.  Two  different  input  deflection  vari¬ 
ables  are  available  for  the  x  and  z  axes.  BZP  and  BXP  should  be  used 
when  absolute  ouputs  (referenced  to  Inertial  spacr  '  are  desired;  DZP  and 
DXP  are  used  when  outputs  relative  to  the  moving  platform  are  desired. 
The  equations  of  motion  are  in  Appendix  C.  Background  on  their  deriva¬ 
tion  and  development  are  given  in  a  series  of  prior  reports  (References 
1-3,  6,  10). 
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USING  BI00TB-80 


Three  job  steps  an  required  for  a  complete  run  of  BIODYN-80.  The 
CREATE  program  is  used  interactively  via  a  remote  terminal,  to  assemble 
the  PARAMETER  and  CHOICES  files.  These  files  are  then  read  by  BIODYN  and 
the  transfer  functions  computed,  stored  on  the  TF  file,  and  listed  at  a 
line  printer.  The  TF  file  is  in  turn  read  by  the  PLOT  program  to  pro¬ 
duce  any  Bode  plots  desired.  Each  of  these  steps  is  described  in 
detail.  Note  that  all  user  Interaction  with  the  computer  assumes  the 
user  is  logged  in  to  CDC'a  INTERCOM.  It  is  assumed  that  the  user  has 
access  to  Intercom  set  "AFMLH  where  the  program  and  its  files  reside. 
For  other  setname  groups,  please  consult  your  Computer  Center  Represen¬ 
tative  to  get  on  the  AFML  list,  or  for  instructions  on  how  to  use  a 
"modified  Job  stream  listing"  (Reference  8). 

A.  CREATE 

Because  of  the  large  number  of  parameters  required  by  BIODYN-80,  it 
was  decided  that  a  dedicated  Interactive  "editor"  for  creating  parameter 
files  was  needed  to  produce  rational,  error-free  input  files,  and  there¬ 
by  avoid  aborting  during  the  subsequent,  and  expensive,  batch  run  of 
BIOOYN  (for  example,  due  to  a  missing  decimal  point).  Furthermore,  this 
approach  allows  the  maintenance  of  a  "catalog”  of  PARAMETER  files  for 
various  postures  and  situations,  from  which  the  most  appropriate  one  can 
be  easily  selected  and  modified  for  a  particular  situation.  The  CREATE 
program  serves  in  this  regard  by  allowing  a  user  to  interactively  assem¬ 
ble  new  parameter  files  or  modify  existing  ones.  CREATE  also  performs 
various  data  checks  and  identifies  common  parameter  entry  errors.  This 
section  details  the  use  of  the  CREATE  program. 

1 .  Entering  CREATE 

CREATE  makes  use  of  five  files,  as  defined  by  the  PROGRAM  card: 

■  INPUT  file  (the  user's  terminal) 
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TAPE4 


TAPE  5 


OUTPUT  file  (the  user's  terminal) 


T 


TAPE7  »  CHOICES  file  (if  new  one  is  created) 

TAPE8  -  PARAMETER  file  (new  or  modified) 

TAPE 20  -  PARAMETER  file  (existing) 

The  PROGRAM  card  automatically  assigns  the  INPUT  and  OUTPUT  files  to  the 
appropriate  tape  units.  If  an  existing  PARAMETER  file  is  to  be  ac¬ 
cessed,  it  must  be  assigned  prior  to  executing  CREATE.  The  job  stream 
to  be  input  by  the  user  looks  like  this  (detailed  examples  are  given 
later,  in  Appendices  D  and  E): 

REQUEST, TAPE 7, *PF 
REQUEST, TAPE8,*PF 
ATTACH, TAPE20,pfn^ 

ATTACH, EXECRT 
EXECRT 


(pfn^  is  the  PARAMETER  file  name 
selected  by  the  user  from  the 
permanent  file  catalog  within 
the  setname  AFML  group) 

I 


Now  the  user  is  at  the  entry  level  of  CREATE.  The  first  issue 
addressed  is  creating  the  PARAMETER  file,  or  modifying  an  existing 
PARAMETER  file.  The  computer's  first  prompt 

NEW  FILE? 

asks  the  user  to  identify  which  option  is  desired.  Answers  are  YES  or 
NO;  any  other  response  will  cause  the  query  to  be  repeated. 

2.  Modifying  a  PARAMETER  File 

If  the  answer  is  NO,  the  values  stored  on  the  file  assigned  to 
TAPE20  are  read.  The  user  may  list  those  values  by  responding  to  the 
next  prompt 

LISTING  DESIRED? 
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with  YES.  A  NO  response  will  skip  over  the  listing  option.  Any  other 
response  will  cause  repeat  of  the  question.  The  listing  produced  can 
either  be  long  or  short.  The  long  listing  provides  parameter  defini¬ 
tion,  mnemonic,  units  and  current  value,  while  the  short  listing  prints 
only  values,  grouped  according  to  each  element  (see  Table  A-2,  Appen¬ 
dix  A).  Examples  of  each  are  given  in  Figures  7  and  8,  respectively. 

Response  to  the  prompt: 

(LO)NG  OR  (SH)ORT  LISTING? 

determines  which  listing  to  generate.  A  response  other  than  LO  or  SH 
will  cause  repeat  of  the  question. 

Values  are  now  ready  for  change,  as  indicated  by  the  message: 

INPUT  MNEMONIC  (5  CHARS)  AND  VALUE  FOR  EACH  CHANGE, 

TO  EXIT,  TYPE  XXX. 


To  change  any  desired  parameter  value,  its  mnemonic  from  the  second 
column  of  the  definitions  in  Appendix  A  is  input.  Mnemonic  must  be  left 
Justified  within  a  5  character  field  followed  by  a  comma,  followed  by 
the  new  value.  The  mnemonic  is  then  checked  against  the  internal  list 
of  permissible  mnemonics;  e.g.,  if  ABDCE  is  not  found,  the  message: 

MNEMONIC  ABCDE  NOT  PERMISSIBLE.  PLEASE  REINPUT. 

appears  and  the  user  should  try  again.  Likely  errors  include  typos,  and 
forgetting  to  pad  the  5  character  field  with  trailing  blanks. 

Once  a  correct  mnemonic  is  identified,  its  new  value  is  compared 
with  the  recommended  range  of  values  for  that  parameter.  If  the  value 
is  not  within  that  range,  a  warning  message  appears: 

WARNING  -  RECOMMENDED  RANGE  FOR  THIS  PARAMETER  IS 
0.5  TO  3.5  GS 

DO  YOU  WANT  TO  CHANGE  THE  VALUE? 
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Figure  7  Example  Long  Format  File  Listing 


20 


■  ’  '  :  K  r1.  V-  • 

XMS 

Kb 

T .  00000 

‘t«:  ■ 

BS 

N/h '  S 

? 4 0.000 

■;.  :  j  !  !■  1.  T"  .if!  i  t  H  f 

XKS 

N/M 

J.  0000 .0 

1  ;  1  ;  I  :  Tl!  i  !  ’.iJI  :  ; 

» iRiVilf  1  h.K* 

CKS 

.... 

I.  .  00000 

Hi  1 1  • 1. 1 i  11  WL.  , t  MS  i  1  i. 

1 l V  10  NX 

sx 

N/G 

0  ♦ 

Bt  -  HUt  ■  ■  i ;  i  SI-  r>  •  >  i  i  1  v 

1 1 f  TO  NZ 

sz 

N  /  b 

0. 

:  ;  :  r.  .  .1  1: 

THE 

DEG 

/O * oooo 

•  i  -  *  i  !•  1-!  .  !  !  ! 

XLC 

M 

. 55 00 00 

LA: Phi  •!  Hi  1 

i  rio  i  uk 

XKSC 

/M 

.  TOoOOOli?  -O' 

r  ■■  Pi:  h  i.ifTHPI  !•'  ■  i- 

OR  MAE > 

BAR 

N/M/S 

0 . 

v  i  ?  1.  !•  {■  ness 

< NOR HAL ) 

XKAR 

N/M 

0 

V-  !■  ,  .  1  : :  1 1  'I  f  !  P  '  T 

ANA  > 

BA  t 

N/M/S 

0  v 

•  V  11  'h-h!  ;  T  i  j-| 

i  f  AN(T ) 

XKAT 

N  /  M 

0  ♦ 

'  i-.  !  :..i  i  }  •.  H  ;  \  f 

tl  !  ; i  (  AN;  i::I 

XL!?  R 

M 

■-.)  * 

'  .  \  }•■'  M  W!  !  •  )  t.M'V' 

1 1 ; ' n  1 : i- I 

AKMR 

.... 

0  ; 

• ! v  !  :  r  } !  !•  OK  i  ;  IaA! 

ANH  1-  A.D. 

NM 

- 

()  v 

VL.  ;  il  .  '.!•  !  I  . !  v  n  S 

:  1.1:  1  .1.  ."HR NTS 

s  :i. 

•• 

GO , 0000 

■1  .  |»  '  '  •  ;  '  •  !f\*  |  ,  f  » 

ICC 

S 

V  *;  i  V  J()  01“  J  j 

■J  :  :  i;  t  r  '  ’  ."i  :  |  :  if 

Oi  l  A  t 

res 

♦  S90000E-0.I  j 

.  •!  .  .  !■ 

X  K  A  A 

N/N 

1  ♦  ooooo 

i  i"!  ,i  !  ;  .1  ,  i  ?!  v 

ZA 

.... 

. SO 00 00  ! 

■1.1  1  i.ir  hi  i!!l  :  J(  ■; 

IN  1-  M,(l. 

W  A 

R/S 

16, 0000  ! 

i  ill-  ill  -  i  1 1  i  IN  1  .  A  , 

0 .. 

I  A  A 

S 

0 .  | 

-1  1  :  :  1  !iUI  ''iiiil'l  I.' 

BM 

N  -Ti/  3 

1  •>  00000 

ii  h  lii'  i  AI  ■>  h  A  r 

EM 

XKM 

N/M 

2 . 00000 

Si:  R  >  1  !'  1  A  -1  r  l.i  1.  i:: 

HEN T  TAMPER 

BE 

N/M/S 

2 . 43000 

'■■:  -  i  ,  V-.  ;  •  1"  I-  |  f 

HI  N  i  OR  A  ft 

XKE 

N/M 

40.0000  ! 

!"•■  i;  J  il  .1  . 

B  r 

N/M-'S 

0.  | 

i  r  :■■■!.'  LI  I  !  ;'■  A  ■ 

XKT 

N  /M 

80,0000  1 

:  "III.  1  Ml  All  i  M.'l  ;  r-! 

xi;sp 

N  /  h 

5.00000 

■•■■■  ,  ■  ;m  :  * i:  ah  i  i 

i  HNS  TAN  I 

ISP 

S 

,9090001?.: -o.i.  ; 

■ :  •  1.1 ! A  !  !  ,.•■■■  i  !  hi" 

iT  NS  T  ANT 

TPS 

C; 

0,  1 

i  ".‘III  !  'll  i  T  !'  1  I- 

1 1  Ti  l.C. 

TSS 

b 

0,  ] 

; . . i  I  .  I  ;  "  i  "  ;  Hi  i  i  M 

XKG 

N/N 

,500000 

."Ml  1  -  ■  III  (  ihl-  S'.h 

v\  S  ?  t)  N  I 

TG 

i:> 

,  550000E  ~0 .1  j 

=  -  ■  !  ''ll  :  "■  i'  1  ’  !  I  111 

1  i  Hi  OONST 

r  z 

s 

\J  * 

■  ;  .  .• .  l  '  hi  ■  hf 

S  1  ANT 

TP 

c;; 

0,  1 

•  •  i  1  '  i 

;  .1 

XKBE 

R  /  R 

7.07950  | 

XKRE 

R/R  /S 

0 . 

i ION  LOOP 

TV 

s 

,  430 000 1?  ■  o I 

•  •  ’  .  :  i  •  •  r 

'•1.1  l:  PEi.i 

AL  PH  A 

R/S 

.300000  ; 

.  ! 

- 1 

XKP 

R/R 

\t  ♦ 

t 

■  :  H  SI  NSR 

XKIE 

.... 

.  670000  ; 

!  ifi  SR  NSR 

XKVli: 

- 

,100000  i 

re? 

5 

,  1 000001?  -0  1.  , 

■  ■  . 

■ '  ; 

TEM 

o 

.100000  1 

■ 

:  n/T  i 

TEN 

s 

,  1230001?  "  0  I  ; 

. .  if*  !  i  i 

1  HOC!' 

OE 

1 .00000 

i  :  ‘  ■  •  :  ,  •  1 

IE 

.... 

, 630000 

•  .  t  •  • 

I  ■  : 

WE 

R/S 

3  i.  .•  .2 0  | 

I 

Figure  7 

(Concluded) 

1 

—>r«  - .  -*  ' 


PARAMETERS  FOR 
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Figure  8  Example  Short  Format  PARAMETER  File  Listing 
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If  a  mistake  has  been  made  In  typing  the  new  value,  a  YES  answer  will 
allow  retyping  of  both  the  mnemonic  and  Its  corrected  value.  If,  how¬ 
ever,  the  value  Is  correct,  even  If  out  of  range,  e  NO  answer  will  cauae 
the  value  to  be  stored  as  It  stands.  Any  response  other  then  YES  or  NO 
will  cause  the  warning  message  to  be  repeated. 

The  range  of  parameters  given  in  Appendix  A  la  ee  wide  as  can  be 
allowed  for  typical  applications  of  BI0DYN-80  (e.g.,  arm  joint  limits). 
Values  outside  this  range  are  used  at  the  user's  risk,  and  precise  pos¬ 
tural  angles  for  the  given  case  should  be  carefully  checked  at  this 
point.  (It  la  good  practice  to  check  your  lengths  and  angles  by  drawing 
a  stick  figure  to  scale,  using  the  desired  values.) 

When  all  the  desired  changes  have  been  made,  an  Input  of  XXX  causes 
control  to  leave  the  modify  mode.  An  opportunity  to  change  the  case 
title  is  then  presented  as 

TITLE  IS  STD  CREWMAN,  HANDS  ON  KNEES 

CHANGE  DESIRED? 

A  YES  answer  causes  the  message 

NEW  TITLE: 


to  appear  and  a  new  title  (60  characters  maximum)  can  be  typed.  A  NO 
response  skips  this  section;  other  responses  cause  the  prompt  to  repeat. 

After  all  desired  changes  have  been  made,  the  user  is  once  again 
given  the  opportunity  to  list  the  file,  using  either  the  LOng  or  SHort 
format,  as  described  above.  At  this  point  the  computer  writes  the  new 
values  to  the  PARAMETER  output  file,  TAPE8. 

3.  Assembling  a  New  PARAMETER  Pile 

If  most  of  the  parameters  are  new,  then  instead  of  modifying  an 
existing  file  the  initial  prompt 
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v  -  jsm-  • .  - 


NEW  FILE? 


Is  answered  YES.  Every  parameter  value  must  then  be  input.  The  message 

INPUT  VALUE  FOR  EACH  MNEMONIC.  FOR  FURTHER  EXPLANATION,  TYPE? 

Initiates  this  process.  For  each  parameter,  the  mnemonic  is  printed  as, 
for  example.  Viewing  Distance,  VD: 


VD  - 

and  the  value  is  accepted  in  floating  point  format.  If  the  user  is 
uncertain  as  to  the  range  of  appropriate  values,  an  expanded  prompt  may 
be  requested  by  typing  ”?”.  The  expanded  prompt  for  VD  is,  for  example: 

VD  -  T_ 

VD  VIEWING  DISTANCE 
RECOMMENDED  RANGE  -  0.5  TO  3.0  M 
VD  - 

the  user  can  then  enter  the  desired  value. 

As  each  value  is  read,  it  is  checked  against  the  recommended  range 
of  values,  in  order  to  weed  out  input  errors  due  to  wrong  sign  or  incor¬ 
rect  decimal  placement.  If  the  value  is  not  within  range,  the  warning 
message  given  previously  appears.  The  user  always  has  the  option  of 
modifying  or  saving  any  value  typed. 

This  process  continues  until  all  ?6  parameters  have  been  entered. 
The  case  title  is  then  added  as  a  response  to 

INPUT  TITLE  FOR  THIS  CASE: 

At  this  point,  a  LOng  or  SHort  format  listing  may  be  generated  as  dis¬ 
cussed  above.  The  file  is  then  written  by  the  computer  to  the  PARAMETER 
file  designated  TAPE8. 


4.  Assembling  the  CHOICES  File 

The  next  section  of  CREATE  is  used  to  assemble  a  new  CHOICES  file 
which  specifies  the  transfer  functions  to  be  output.  Note  that  there  Is 
no  provision  for  modifying  an  existing  CHOICES  file;  any  changes  to  an 
existing  CHOICES  file  must  be  made  by  creating  an  entirely  new  file. 
The  prompt 

NEW  CHOICES  FILE? 


has  two  responses:  a  NO  allows  the  user  to  skip  this  entire  section  If 
an  existing  CHOICES  file  will  be  used  as  Input  to  BIODYN;  a  YES  response 
will  cause  prompting  for  the  CHOICES  file  components  as  described  below. 
Any  other  response  will  cause  the  prompt  to  be  repeated. 

A  YES  response  to  the  next  prompt 
BIODYN  TFS  DESIRED  FOR  PIVIB? 


will  automatlcelly  generate  the  ten  transfer  functions  to  be  Included  In 
the  PIVIB  Input  file.  This  Is  mandatory  if  Interaction  with  PIVIB  is 
anticipated.  The  response  and  forcing  function  variables  Involved  are 
shown  later,  in  Figure  12.  A  NO  answer  skips  this  section;  other 
answers  cause  repeat  of  the  question. 

All  other  desired  transfer  functions  are  entered  by  the  user  in  re¬ 
sponse  to  the  following  computer-specified  format: 


TRANSFER  FUNCTION  INPUT: 

FIRST  LINE  -  RESPONSE  MNEMONIC,  FORCING  FUNCTION  MNEMONIC 
(AAA,  AAA);  ENTER  XXX  TO  STOP 
SECOND  LINE  -  PLOTTING  INFORMATION,  5  ITEMS: 

BODE  LOWER  FREQ.  LIMIT 

BODE  UPPER  FREQ.  LIMIT 

BODE  UPPER  PHASE  LIMIT  (0.  DEFAULTS  TO  200.) 

BODE  LOWER  PHASE  LIMIT  (0.  DEFAULTS  TO  -400.) 

LIST  (1.  TO  LIST  TABLE,  0.  FOR  NO  LIST) 

IF  NO  PLOT  DESIRED,  ENTER  0.  FOR  ALL  ITEMS 
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Note  that  the  parameters  for  each  transfer  function  are  entered  on  two 
consecutive  lines,  each  Item  separated  by  a  comma,  each  line  followed  by 
CR. 


The  mnemonics  requested  are  three-character  codes  as  defined  In 
Appendix  B.  The  program  checks  each  mnemonic  against  the  list  of 
acceptable  codes.  If  a  match  is  not  made,  the  error  message 

XYZ  NOT  PERMISSIBLE,  PLEASE  REINPUT, 
is  given  and  the  entire  first  line  must  be  retyped. 

If  Bode  plots  are  desired,  the  frequency  and  phase  limits  must  be 
entered;  if  not,  zeros  are  assumed.  The  frequency  limits  must  bound  a 
frequency  range  of  no  more  than  3  decades.  If  frequencies  with  an 
unacceptable  range  are  entered,  the  message 

MAX  FREQUENCY  RANGE  IS  3  DECADES 

PLEASE  REINPUT  ENTIRE  LINE 


appears  and,  all  five  items  for  that  transfer  function  must  be  re-typed. 
The  phase  limits  are  unbounded,  but  entering  zeros  will  cause  default  to 
+200.0  deg  and  -400.0  deg  for  the  upper  and  lower  limits,  respectively. 
If  a  tabular  listing,  including  frequency,  amplitude,  and  phase  at  20 
increments /decade  intervals  is  desired,  a  1.0  should  be  entered  as  the 
fifth  item;  otherwise  enter  0.0. 

The  CREATE  program  stops  when  an  XXX  is  read  in  the  CHOICES  file 
section.  The  two  files  created  must  be  saved  as  permanent  files  by 
giving  the  following  INTERCOM  commands: 

CATALOG ,TAPE7,pfn2,RP-10 
CATALOG, TAPE8 , pfn3,RP-10 

where  pfn2  is  a  new  (seven  character)  permanent  file  name  assigned  to 
the  CHOICES  file  and  pfn3  is  a  new  (seven  character)  permanent  file  name 
assigned  to  the  PARAMETER  file.  The  user  is  now  ready  to  exercise 
BI0DYN. 

The  user  can  interactively  create  a  number  of  PARAMETER  and  CHOICES 
files,  each  uniquely  named  for  subsequent  use  in  the  BI0DYN  program.  In 
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most  cases  one  same  set  of  transfer  function  choices  will  be  adequate 
for  a  variety  of  different  parameter  sets. 

B.  Bloom 

Because  of  memory  size  limitations  presently  Imposed  by  the  INTERCOM 
operating  system,  BIODYN  must  be  run  In  the  hatch  mode.  A  batch  Job  can 
be  submitted  through  INTERCOM  using  the  following  command  sequence 


EDITOR 
CREATE, S 

ABC, CM1 50000, STCSA. 

ATTACH, TAPE 7, pfn2«  User  supplies  desired 

ATTACH, TAPE8,pf no.  file  name  for  each  pfn 

REQUEST, TAPE 19, *FF. 

ATTACH, EXEBIO. 

EXEBIO(TAPE4 .OUTPUT ,TAPE7,TAPE8,TAPE19 .TAPE2I ) . 
CATALOG , TAPE 19 , pf n4 , RP-999 . 

*E0R 

SAVE.GOBIO .NOSEQ 
END 

COMMAND  -  BATCH, GOBIO, INPUT, HERE 


The  first  section  of  this  batch  request  assigns  the  PARAMETER  and 
CHOICES  files  produced  by  CREATE  to  the  tape  units  which  BIODYN  reads. 
The  BIODYN  program  Is  then  retrieved  and  executed- 

BIODYN  first  reads  the  PARAMETER  file  and  uses  It  to  set  up  the 
equations  In  the  form 


Ax  -  Bu 


where  x  is  a  vector  of  response  variables  (currently  dimensioned  48)  and 
u  is  a  vector  of  forcing  function  variables  (dimensioned  9).  The  de- 
sired  transfer  functions  (as  listed  In  the  CHOICES  file)  are  then  evalu¬ 
ated,  using  Cramer's  rule  and  several  advanced  factorization  algorithms 
In  order  to  obtain  the  first-  and  second-order  poles  and  zeros.  These 
are  all  written  to  a  file  called  TAPE 19  (the  TF  file),  as  well  as  to  the 
line  printer  (via  the  OUTPUT  file).  The  format  used  by  the  program  for 
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'tsm+d  m  mm wm-ae*  •- 


printing  the  transfer  functions  is  shown  in  Figure  9,  while  the  foraat 
for  storing  these  transfer  functions  on  TAPE19  is  given  in  Figure  10. 

The  TF  file  (TAPE19)  is  used  by  the  PLOT  program  for  generating  the 
quick-look  Bode  plots  on  the  line  printer.  Thus,  the  last  lines  sup¬ 
plied  by  the  user  in  the  batch  job  stream  "rewind”  this  file  and  assign 
it  a  permanent  file  name  pfn4. 

c.  not 

The  Bode  plots,  if  desired,  may  be  generated  at  the  user's  terminal 
or  written  to  a  local  file  for  later  examination,  as  the  user  wishes. 
The  user's  job  stream  looks  like  the  following: 

ATTACH .TAPE 19, pfn4.  pfn4  is  the  file  name  of 

ATTACH, EXEPLT.  the  plot  instructions 

EXEPLT ( OUTPUT , TAPE 1 9 ) 

The  output  from  this  routine  consists  of  four  items.  First,  the 
input  file  mnemonics  are  printed,  so  the  user  can  Identify  which  plots 
are  forthcoming.  Then,  for  each  transfer  function,  three  items  are 
printed.  The  transfer  function  dipoles  (first  and  second  order)  which 
have  been  cancelled  are  first,  followed  by  the  numerator  and  denominator 
of  the  transfer  function.  Then,  the  Bode  plot  is  "drawn."  Finally,  if 
a  listing  was  requested  in  the  CHOICES  file,  it  follows  the  Bode  plot. 
An  annotated  example  is  given  in  Figures  11a,  b,  and  c. 

D.  UrmPAOE  WITH  PITH 

Ten  of  the  possible  BIODYN-80  transfer  functions  can  be  used  as 
input  to  PIVIB.  These  are  listed  in  Figure  12.  The  input  file 
structure  for  PIVIB  is  quite  lengthy  and  complex,  involving  over  200 
parameters.  An  explanation  of  it  is  beyond  the  scope  of  this  user's 
manual,  and  the  interested  reader  is  referred  to  Reference  7,  the  PIVIB 
User's  Manual.  A  brief  outline  of  the  steps  used  to  generate  the  PIVIB 
input  file  is  given  below. 
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Each  transfer  function  is  to  b»*  interpreted  as 


^output 

Transfer  Function  =  VriEut -  (*) 
A 


‘H5  »£)}•{■»(  «?*  )} 

|lt(*  *  ?;)  •{"(■’*  Vp'^)} 


where  l  denotes  a  product  of  first-  or  second-order  roots,  and  *  Is  the  so-called  "root  locus  gain"  (of  the  hlgh- 
frequency  asymptote): 


High  Frequency  Gain  of  Numerator 
High  Frequency  Gain  of  Denominator 


and  the  various  first-  and  second-order  poles  and  zeros  are  indicated  below.  The  transfer  function  can  also  be 
interpreted  in  Bode  format,  as: 


Transfer  Function 


’)Hn(|*y  *  ')\ 

{n(V*  i)}.|n(^*2£.  i)j 


where  is  the  gain  of  the  low-frequency  asymptote  ("Bode  gain"): 


N)l 


Low  Frequency  Gain  of  Numerator 
Low  Frequency  Gain  of  Denominator 


1 1 -MAR-80  17106 


Date,  Time 


CASE : HEAD  » NECK  »  TORSO  ONLY 


Case  Identifier 


. 10464E-12 

(  46.334  )  (  100.00  )  (  530. A3  > 


<  •: 

.40231 

8.3929 

.  3.7123 

.  7.5273 

>> 

<  ( 

.42607 

10.483 

»  4.4AA3 

#  9.4335 

)  ) 

( ( 

.13834 

17.112 

*  2.199A 

9  16.970 

)  ) 

( ( 

. 97453F -0 1 

32. 74A 

*  3.1912 

*  32.590 

)  ) 

( ( 

. 1 138A 

44.530 

»  5.0704 

f  44.241 

)  ) 

<  < 

.65021 

31  A. 83 

»  20A.02 

,  240.70 

>> 

.  I  2459E4- 1  2 


TRANSFER  FUNCTION  COMPONENTS: 

- High-frequency  gain 

(denominator)  Kp 

-  First-order  poles  ( 1 /T) 


Second-order  poles 

(U,  a,  £id,  u>yr t" 


Low-frequency  gain  (denominator) 


numerator:  KH0/r'7F‘' 


( * 

,00000 

\  . 00000 

)  <  46 

384 

)  ( 

100. 

( 

<  < 

534.5? 

.44231 

*  8.3929 

*  3.7123 

F 

7.5273 

>  ) 

(  ( 

. A  7357 

.  17.729 

f  11.94? 

* 

13.104 

)  ) 

<  ( 

-.35676 

.  35.405 

. -12.631 

F 

33.075 

)  ) 

<  < 

. 29926 

.  44.192 

»  13.223 

F 

42,166 

)  ) 

( ( 

.  A5024 

.  316.83 

f  206.02 

F 

240.70 

)  ) 

- 

.  1 1284EM  2 

-  •  905A5 

- - 

K 


Numerator:  output /input 

<High-freq.  galn>  numerator  alone  Kgj 

<<High-freq.  galn>>  for  transfer  function 


First-order  zeros 


K  -  Kg/Kp 


Second-order  zeros 

«Low-freq.  gain,  transfer  fUnctlon>> 


V 


<Low-fTeq.  gain,  numerator  alone > 


Figure  9 


Format  for  BIODYN-80  Transfer 


Function  Printouts 
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-  .'I'lwiwiim'  win  ii.iai .  n 


Format 


Comment 


Case  title  (15A4) 
D$$ 


NP1,  NP2,  Kp,  0.,  0., 


1/T 

1/T 


Pi 

P2 


o.. 


o.. 


0. 


(Specifies  denominator) 

(214,  G14.6,  5F7.1) 

First-order  poles  (5X,  G14.6) 


4P1,  “pi  Second-order  poles  (ci  >  t*>i) 

CP2,  o>p2  (5X,  2F14.6) 


xxx/xxx 

NZ1 ,  NZ2,  Kz,  WI,  MF,  XU,  XL,  XP 

1/TZ1 

i/tz2 


4*1 1  '°*1 
4*2*  “*2 


XXX/XXX 


First  transfer  function 
numerator 

(214,  G14.6,  5F7.1) 
First-order  zeros  (5X,  G14.6) 


Second-order  zeros  (c2»  “2) 
(5X,  2G14.6) 

Second  transfer  function 
numerator 


where 


NPl  “  No.  first-order  poles 
NP2  ■  No.  second-order  poles 
Kp  -  High-frequency  denominator 
gain 

NZ1  •  No.  first-order  zeros 
NZ2  ■  No.  second-order  zeros 


K  -  High-frequency  numerator  gain 
WI  ■  Bode  plot  lower  frequency  limit 

WF  ■  Bode  plot  upper  frequency  limit 

XU  ■  Bode  plot  upper  phase  limit 

XL  ■  Bode  plot  lower  phase  limit 

XP  ”  Listing  switch 


Figure  10,  Format  Used  for  BI0DYN-80  TAPE 19  File  of  Transfer  Functions 
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BIODYN  T.F.  05/20/80 

SEATED  CREWMAN  -  VIEWING  DISTANCE 


Title 


15.13.57. 
--  9  CM 


0«* 

RHD/DZP 


Dump  of  f ile 


DIODYN  T.E. 


BIODYN  T.F.  05/20/00 

SEATED  CREWMAN  -  VIEWING  DISTANCE 


15.13.57. 
9  CM 


Title 


RHD/DZP  - 

D** 

FIRST-ORDER  DIPOLES  CANCELLED 
46.  384000 
100.00000 

SECOND -ORDER  DIPOLES  CANCELLED 
.44230900  .  8.3929300 

.40022000  .  50.002500 

.65024400  »  316.83100 


First  transfer 
function 

Cancelled  dipoles 
(only  exact  can¬ 
cellation 


numerator: 


-.81361 


(0 

» 

)  (0. 

)  (  10 

.955 

>  ( 

11.950 

< 

22.500 

)  (  107.75 

)  <  53 

5.79 

) 

<  < 

.22639 

.  a. 4561 

f 

1.9144 

f 

8.2365 

>>  ^ 

<  < 

. 78552 

»  13.624 

r 

10.702 

f 

8.4311 

>> 

<  ( 

.82741 

.  20.272 

f 

16.777 

T 

11.379 

>)  > 

<  ( 

-.16669 

f  21.651 

-3.6090 

f 

21.348 

>>  ( 

<  < 

♦  2273? 

.  47.017 

y 

10.700 

f 

45.783 

>>  J 

DENOMINATOR 

: 

( 

10.957 

)  (  11.983 

)  <  22 

.499 

)  ( 

107.81 

( 

531.02 

) 

(  / 

. 22869 

»  8.4602 

t 

1.9366 

r 

8.2437 

(  ( 

.40550 

r  13.701 

f 

3  ♦  5538 

r 

12.524 

>>  ) 

(  ( 

.24026 

>  16.682 

r 

4.0081 

r 

16.193 

>>  ( 

(  ( 

.82362 

f  20.305 

r 

16.723 

f 

11.516 

>)  \ 

(  < 

•94299F- 

01.  31.293 

t 

2.9509 

f 

31 . 153 

>> 

(  ( 

. 19289 

.  46.973 

r 

9.0608 

9 

46  091 

)>  ’ 

■  .  11057F- 

02  - 

— 

>h- 


High  frequency  gain 
First-order  zeros 

Second-order  zeros 


First-order  poles 
Second-order  poles 


Low  frequency  gain 


r:  1  n  M  I  N  I  .  F  . 


MOrifN  r  ,F .  05/20/80  15*13*57* 

->i.ah;d  crewman  -  viewing  distance  *  ?  cm  *- 


RHD/DZP 

n$* 


XXX  MAGNITUDE  (DEO  - 

20.  0.  -20.  -40. 


-60. 


I -'Hi. 


Vj. 


x  Magnitude 


-80. 


100. 


Title 

Top  scale 

XXX  is 
magnitude 


29. 


20. 


].  4, 


I  0. 


I 


5.0 
4 . 0 

2  >*i. 


i  .  1 


i  ■  0 

}<■/  s 


X  0 

ox 

0  X 

X 
X 

X 

X 


1 

Phase 

I 


o 

x  □ 
x  o 
x 

ox 

0  x 

0  X 
0  x 
0 
0 
0 
o 
0 
0 
0 
0 
0 
0 
0 
0 
0 


+360°  phase  shift  for 
o  for  points  below  -l8o° 


Values  at  ">"  (in  rad/sec) 
and  at  — -  (in  Hertz)  are 
exact . 


1.00. 


0. 


X 

-100.  -£00.  -300.  -400. 

000  =  PHASE  <HEG>  - 


Bottom  scale 
for  000  is 
phase 


b)  Printer -Genera ted  Frequency  Response  (Bode  Plot) 
Figure  11  (Continued) 


0 
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Title 


Listing  of  frequency 
(rad/sec),  amplitude  (dB) 
and  phase  (deg)  at  20 
evenly- spaced  increments/ 
decade.  (This  listing  is 
only  generated  if  XL  *  1.0 
in  the  CHOICES  file  for  this 
transfer  function) 


c)  Listing  of  Frequency  Response  Values 


Figure  11  (Concluded) 
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Note:  i>  denotes  a  typed  blank  space 


Output 

BIODYN  TF 

PIVIB 

TF 

Shoulder: 

JiXS 

SB 

Chi 

Inertial  horizontal  shoulder  deflection 

(S) 

BXP 

Chi 

S 

Horizontal  platform  deflection 

bxs 

SB 

Chi 

Inertial  horizontal  shoulder  deflection 

BZP 

Ch2 

s 

Vertical  platform  deflection 

tzi 

■i 

Ch3 

Inertial  vertical  shoulder  deflection 

BXP 

Chi 

s 

Horizontal  platform  deflection 

JSZ1 

Ch3 

Inertial  vertical  shoulder  deflection 

BZP 

CH2 

s 

Vertical  platform  deflection 

Head: 

DTH 

a 

Ch3 

Inertial  head  pitch 

(H) 

BXP 

Chi 

H 

Horizontal  platform  deflection 

DTH 

m 

Ch3 

Inertial  head  pitch 

BZP 

Ch2 

H 

Vertical  platform  deflection 

DZH 

m 

Ch2l 

Inertial  vertical  head  deflection 

BXP 

Chi  1 

H 

Horizontal  platform  deflection 

DZH 

m 

Ch2 

L 

Inertial  vertical  head  deflection 

BZP 

Ch2 

Vertical  platform  deflection 

Stick 

let 

m 

Chi 

Longitudinal  stick  deflection 

(C) 

DXP 

Chi 

lc 

Horizontal  platform  deflection 

tet 

m 

Chi 

Longitudinal  stick  deflection 

DZP 

Ch2 

c 

Vertical  platform  deflection 

Assumptions:  Of  the  3  vibration  Inputs  possible,  2 
are  X  and  Z.  Assume  X  is  Channel  1,  Z  is  Channel  2 
for  the  above  table. 


Figure  12  BIODYN-80  Transfer  Functions  Which  May  Be 
Used  as  PIVIB  Input 


1.  Specify  title  for  run 

2.  BDMOD  -  enter  the  blodynaolc  response  module 

a.  FREQS  -  specify  all  frequencies  required  for 
blodynamic  analysis 

b.  VIBR  -  specify  vibration  input/source  transfer 
functions 

c.  BIOTR  -  specify  blodynamic  transfer  functions. 
This  Is  the  Interface  between  BI0DYN-80  and 
PIVIB. 

d.  BDOUT  -  specify  blodynamic  outputs  for  printing 

e.  BCOMP  -  given  all  Information  specified  above, 
compute  and  print  blodynamic  quantities 

3.  PVMOD  -  enter  the  pllot/vehicle  module 

a.  Specify  various  P/V  parameters  In  state  vector 
form 

b.  VBINT  -  controls  communication  between  BDMOD  and 
PVMOD 

C.  PCOMP  -  perform  pilot/vehicle  computations  and 
print  results 

d.  FDREP  -  compute  and  print  frequency  domain  mea¬ 
sures,  e.g.,  describing  functions,  remnants, 
signal  spectra 

4.  Repeat  2  and  3  for  each  case 

5.  END 


B.  FOB  ASSISTANCE 


New  and  inexperienced  users  should  have  an  experienced  INTERCOM  or 
batch  user  assist  them  throught  the  first  application.  For  technical 
questions  concerning  this  program  or  models,  please  call  or  write: 
Henry  R.  Jex  at  Systems  Technology,  Inc.,  13766  S.  Hawthorne  Blvd., 
Hawthorne,  CA  90250.  Telephone  Number  (213)  679-2281. 
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APPENDIX  A 

DEFINITION  OF  PARAMETERS,  SYMBOLS,  RANGES  AND  TYPICAL  VALUES 

TABLE  1  PARAMETER  AND  SYMBOL  DEFINITIONS,  RANGES,  AND  TYPICAL  VALUES 


PARA¬ 

METER 

FORTRAN 

tm 

UNITS 

LOCATION 

DEFINITION 

SOURCE* 

(Refs.) 

- REASONABLE - 

- 1 

TYPICAL  VALUES  (STD  CASES)  | 

RANGE 

MINIMUM 

MWIMJM 

SlDPLTt 

t 

STDTCRW 

LOCR  BODY 

* 

m 

kg 

Lower  body  ness  (Hips) 

T(l,2,3) 

10. 

25. 

14. 

14. 

Bb 

B8 

nAn/s 

Lower  body/seat  cush- 

F(l,4,5) 

100. 

2,000. 

1290.1 

315. 

ion  damper 

XKB 

NAn 

Lower  body/seat  cush- 

F(l,4,5) 

20,000. 

100,000. 

29585. 

49360. 

ion  stiffhess 

KfiS 

aa 

NAn 

Series  spring  g-adi- 

F(l,4,5) 

50,000. 

200,000. 

71519. 

119322. 

ent  in  lower  body 

0L 

TH. 

deg 

Seat  (hips)  tilt 

M(l) 

-20. 

+80. 

13. 

13. 

angle 

G 

G 

m/s2 

Trim  vertical  accel- 

0.1 

100. 

9.8 

9.8 

a- at ion 

TORSO 

Mr 

XMT 

kg 

Torso  mass 

T(l,2,3) 

10. 

30. 

18. 

18. 

it 

ZIT 

kg.ntf 

Torso  inertia 

T(  1,2,3) 

0.5 

2. 

0.8 

0.8 

&TB 

BFB 

N.m 

Torso  damping 

F(l,4) 

10. 

100. 

16.75 

40. 

KTB 

XKTB 

r®/5 

N.m 

Torso  stiffness 

G,F(1,4) 

100. 

2,000. 

500. 

500. 

* 

r3 

9t 

THT 

*g 

Angle  of  XLT  fron 

M(  1,3,4) 

-30. 

+90. 

10. 

-15. 

vertical 

ILt 

XLT 

m 

Torso  c.g.  to  hips 

T,M(1,2,3) 

0.12 

0.40 

0.15 

0.15 

pivot  length 

Pin 

THTN 

deg 

Angle  of  XLTN 

M,T{1,2,3) 

-30. 

+90. 

5. 

-5. 

|lTN 

XLTN 

m 

Torso  c.g.  to  neck/ 

T,M(1,2,3) 

0.15 

0.30 

0.3 

0.3 

pivot 

P  s 

1HS 

deg 

Angle  of  XLS 

M,T(1,3) 

-30. 

+90. 

5. 

-5. 

H-s 

XLS 

m 

Torso  c.g.  to  shoul- 

0.10 

0.30 

0.3 

0.3 

der  length 

f’SIDPLl 

r'  =  STD  pi  1 

ot,  stiff  centerstick 

+"SDTCW  =  SID  craman,  hand  on  knees 


♦Sources  Le^nd:  Measured  (e.g.,  photo)  or  specified  by  situation  (continued) 

Tabulated,  cadaver  or  anthropometric  tables 
Fitted,  to  transfer  function  data 
Ties  st  irate  (on  basis  of  physical  properties) 
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TABLE  1.  (Continued) 


- ES50TOBK - 

1 - 

1 

PARA- 

FORTRAN 

UNITS 

LOCATION 

SOURCE* 

RANGE 

TYPICAL  VALLES  (STD  CASES)  | 

METER 

name 

CEFINITICN 

(Refs.) 

MINIMUM 

MAXIM*!, 

STDPLT 

S1DCRW 

FEO( 

m 

kg 

Neck  mass 

T(l,3,6) 

0. 

3.0 

0.0** 

1.24 

In 

m 

kg. 

Neck  inertia 

T{1,2,3,6) 

0. 

0.010 

0.0** 

0.0024 

Bn 

BNT 

N.m 

r37s 

N.m 

Neck/torso  danper 

G,F(1,5) 

0. 

0.40 

0.25 

0.1 

% 

XKNT 

Neck/torso  stiffness 

G(l,5) 

10. 

80. 

50.0 

50.0 

F3 

9n 

THN 

deg 

Angle  of  XLN 

*1.3) 

-90. 

+90. 

-20. 

-30. 

Ln 

XLN 

m 

Neck  length 

M,T(1,2,3,6) 

0.05 

0.25 

0.1 

.12444 

lN1 

XLN1 

m 

Neck  c.g  to  neck/ 
torso  pivot 

T,G(1,3) 

0.02 

0.15 

0.05 

.06222 

Hi 

XMi 

Hea  mass 

T,G(J,2,3,5) 

2.0 

10.0 

helmeted 

4.34 

3.1 

Ih 

ZIH 

kg  .m2 

Head  inertia 

T,G(1,2,3,6) 

0.010 

0.100 

helneted 

0.039 

0.0303 

Hti 

Ml 

N.m 

rS7s 

N.m 

Head/neck  danper 

G,F(1,8) 

0. 

0.40 

0.0** 

0.126 

Kffi 

XKifi 

Head/neck  stiffness 

G,F(1,4,8) 

10. 

30. 

10.** 

15. 

Fa 

Qi 

Oi 

- 

Head/neck  oonpliance 

*1,8) 

0.0 

1. 

0.0** 

1. 

pa-ameter  (CH  =  0. 
locks  head  on  neck) 

locked 

in  locked 

'Hi 

1 

TH 

deg 

Aigle  of  U.G(p1vot-CG) 

M,T(1,2,3,6) 

-30. 

+150. 

70. 

60. 

ILt 

1 

l9v 

XLH 

m 

Head  c.g.  to  head/ 
neck  pivot  length 

T,M(1,2,3,6) 

0.0 

0.100 

0.0* 

0.025 

TEN 

deg 

Angle  of  line  of 

Sl(flt 

M,T(1,3) 

-90. 

+90. 

-30. 

-30.  | 

1 

1 

IVd 

1 

VD 

m 

Viewing  distance 
fron  head/neck 

*1) 

0. 

1000. 

0.75 

0.6863  I 

1 

1 

1 

(pivot  | 

_ l 

•Sources  Legend:  **Head/Neck  Locked 

Measured  (eTg.,  photo)  or  specified  by  situation 
Tabulated,  cadaver  or  anthropometric  tables 
"Fitted,  to  transfer  firction  data 
TLesst  irate  (on  basis  of  physical  properties) 


V 


,  ■  ***  *  vv  -  -*-* 


TABLE  1.  (Continued) 


7YP10L  WUK  (STD  CASES)  | 

PARA- 

FORTRAN 

UNITS 

LOCATION 

S01FCE* 

RANGE 

PtlfcK 

IW€ 

DEFINITION 

(Refs.) 

mm 

5TCTCT 

rsrcw 

ARMS  (Upper,  Lover) 

mi 

XM1 

kg 

Upper  arm  mss 

T,M(1,2,3) 

1.0 

4.0 

1.372 

1.372 

Hi 

n 

kg.n£ 

Upper  arm  inertia 

T(  1,2,3) 

0.01 

0.05 

0.012 

0.012 

k 

TI 

cfeg 

Upper  arm  angle 

M,T( 1,3,8) 

-30.0 

+90.0 

40. 

15. 

14 

D1 

m 

Upper  arm  length 

( =  4  ♦  4) 

M,T( 1,3,8) 

0.2 

0.5 

0.29 

0.29 

ILjAi 

R1 

mM 

Upper  arm  c.g. 
(factional  distance 

T,G( 1,2,3) 

0.2 

0.6 

0.44 

0.44 

1* 

W 

kg 

Lower  arm  mass 

T,M(1,2,3) 

0.5 

3.0 

1.017 

1.017 

H2 

72 

kgjn2 

Lower  arm  inertia 

T,(l,2,3) 

0.01 

0.10 

0.0152 

0.0152 

k 

7E 

deg 

Elbow  angle 

M,T(1,3,8) 

60.0 

170.0 

145. 

115. 

|4 

DP 

m 

Lower  arm  length 

(  =  4  +  4) 

M,T(1,3,8) 

0.2 

0.5 

0.305 

0.306 

14/4 

R2 

m/m 

Lower  arm  c.g. 
(fractional  distance) 

T,G(1,2,3) 

0.2 

0.6 

0.5 

0.5 

|i9u 

GRIP  INTERFACE 

TIJ 

*g 

Grip  interface  angle 

W) 

-45.0 

+45.0 

0.0 

0.0 

IBi/Kl 

TIF 

s 

Grip  interface  time 
constant 

G,F(1,8) 

0.0 

0.5 

0.01 

0.1 

lcI 

a 

m7N 

Grip  interface  com- 

H,G,F(1,8) 

0.0 

1000.0 

.0000558 

1.0 

1 

pliance  Ki-1 

fixed 

"free" 

♦Soirees  Legand:  Measired  (e.g.,  photo)  or  specified  by  situation 
Tabulated,  cadaver  or  anthrcpometric  tables 
Fitted,  to  transfer  find  ion  data 
Ties st irate  (on  basis  of  physical  properties) 
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TAfiLE  1.  (Continued) 


— 

— 

1 

RANGE 

TYPICAL  VALUES  (STD  CASES)  | 

PPM- 

FCRTRAN 

UNITS 

LOCATION 

SOURCE* 

1 

l€7ER 

tm 

CEFINITION 

(Refs.) 

MINIHJM 

MAXIMJM 

STDPLT 

jwirHM 

STICK 

— 

WS 

kg 

StlckAiaid  mass 
(referred  to  grip) 

M.F(1,8) 

0.0 

1.0 

1.31 

4.0 

Bs 

as 

N/m/s 

Stick  danper 

M,F(1,8) 

0.0 

5.0 

2. 

240.0 

Ks 

XKS 

N/m 

Stick  cradient 

*1,8) 

0.01 

20,000.0 

13000. 

10000.0 

Qcs 

ocs 

- 

Stick  conpllance  para- 

*U 

0.0 

1.0 

1.0 

1.0 

rater  (QCS-0  means 
that  XKS-** 

"flrad" 

•fi*d" 

Sx 

sx 

N/g 

flobveiebt  sensitivity 
to  N* 

F,M,S(1,8) 

-5.0 

+10.0 

0.0 

0.0 

Sz 

sz 

N/g 

Bobweigit  sensitivity 
to  Nz 

F,M,G(1,8) 

-10.0 

+10.0 

0.0 

0.0 

\ 

THC 

(feg 

Stick  aigle 

M(l) 

-20.0 

+120.0 

90. 

90. 

Lc 

XLC 

m 

Stick  length  (pivot 
to  center  of  g-1p) 

*1) 

0.0 

2.0 

0.61 

0.55 

*sc 

XKSC 

Units 

Stick  output  scale 

F,G(1,8) 

0.01 

20,000.0 

13000. (Nfa) 

0.01 

m 

factor  (rescales  stick 
displacement  to  "spec- 

ified"  rnlts  (e.g., 
deg,  cm,  etc.) 

BAR 

BAR 

N/m/s 

Ann  rest  dbnper  (nor- 
ral) 

Ann  rest  stiffness 

H.G(1) 

0.0 

1.0 

o.ot 

t 

o.ot 

t 

kAR 

XK PR 

N/m 

M,G(1) 

0.0 

10,000.0 

o.ot 

O.Ot 

(normal) 

t 

t 

Bf 

BAT 

N/m/s 

Arm  rest  danpa* 
(tangential) 

G,M(1) 

0.0 

150.0 

0.0J 

0.0J 

Kf 

XKAT 

N/m 

Ann  rest  stiffness 

G,H(1) 

0.0 

10,000.0 

o.ot 

O.Ot 

(tangential) 

t 

T 

lBR 

XLER 

m 

Elbow  to  am  rest 
distance 

M(l) 

0.0 

0.5 

O.OJ 

0.0J 

An»R 

pm 

- 

Fraction  of  arm 

G(l) 

0.0 

1.0 

o.ot 

0.C* 

weitfrt  on  arm  rest 

t 

♦Sarces  Legend:  Mwsired  (e.g.,  photo)  or  specified  by  situation 
Tabulated,  cadaver  or  anthropometric  tables 
fitted,  to  trwsfer  f met  Ion  data 
ILesst  Irate  (on  basis  of  physical  properties) 
purest  disabled. 
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TABLE  1.  (Continued) 


PARA- 

l^TER 

FORTRAN 

NM€ 

UNITS 

LOCATION 

DEFINITION 

- 1 - 

SOURCE*  |  REASONABLE 
(Refs.)  |  RANGE 

TYPICAL  VALUES  (STD  CASES)  | 

STD CRT 

NBJROKJSCULAR  SYSTEM: 

Actuation  Elements: 

I 

W 

AM 

_ 

Switch  for  central  or 

1 

M  |  0.0 

1.0 

0.0 

0.0 

force  activation  dy- 

|  arm 

head 

names  0  for  arm  oon- 

trol  1  for  head/neck 

1 

control 

1 

Si 

SI 

- 

Overall  scale  factor 

F,G(1,8)  |  0.0 

200.0 

104.4 

50.0 

Kk 

XK1C 

WN 

Neu-onuscular  actua- 

F,G(1,8)  |  0.0 

1.0 

.00044 

.00 

tion  game 

1 

Toe 

TOC 

s 

Neuromuscular  actua- 

F,G(1,8)  |  0.0 

0.2 

.0909 

.0909 

tion  game 

i 

cs 

TCS 

s 

Neu'onuscular  actua- 

F,G(1,8)  |  0.0 

0.3 

0.089 

.089 

tion  lag 

1 

Ka 

XKAA 

N'N 

Gain  of  force  activ- 

1 

F,G(1,8)  |  0.0 

10.0 

1. 

1. 

ation  dynamics 

1 

^a 

ZA 

- 

Oanping  in  force 

F,G(1,8)  |  0.0 

1.0 

0.8 

0.8 

activation  dynames 

1 

wa 

VA 

rad/s 

Natural  frequency  in 

F,G(1,8)  |  10.0 

20.0 

16. 

16. 

force  activation 

1 

dynamics 

T 

a 

TAA 

s 

Tims  delay  in  force 

F,G(1,8)  |  0.0 

0.10 

0.0 

0.0 

activation  dynamics 

|  • 

I 

_ 

^Sources  Legend:  Meastred  (e.g.,  photo)  or  specified  by  situation 
“Tabulated,  cadaver  or  anthropometric  tables 
"Fitted,  to  transfer  function  data 
Tauesst  irate  (on  basis  of  physical  properties) 
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TAELE  1.  (Continued) 


PARA- 
Ft  IER 

FORTRAN 

NWE 

UNITS 

LOCATION 

DEFINITION 

— marnot — 

RANGE 

SOURCE* 

(Refs.) 

Mm 

mill 

smT 

5TDCFW 

NBJROMJSCULAR  SYSTEM; 
Actuation  Elements: 
continued 

EH 

BM 

NAn/s 

"Hills  law"  danper 

G,F(1,7,8) 

0.0 

20.0 

l. 

1.0 

km 

XKM 

N.m 

Spring  in  neroms- 
cular  system 

G,F(1,7,8) 

0.1 

20.0 

2. 

2. 

Be 

I 

NAn/s 

Series  elastic  ele- 
nent  danper 

F(l,7,8) 

0.0 

10.0 

2.431 

2.43 

Kg 

XKE 

N/m 

Series  elastic  ele¬ 
ment  g-adierrt 

G,F(1,7,8) 

0.1 

100.0 

40. 

40. 

Bf 

BT 

NAn/s 

Tendon  danper 

6,F(1,7,8) 

0.0 

10.0 

0. 

0. 

kt 

XKT 

N/m 

Tendon  g-adient 

G,F(1,7,8) 

0.1 

300.0 

80. 

80. 

Ksp 

XKS> 

N/m 

Hiscle  spindle 
model  gain 

G(l,7,8) 

1.0 

10.0 

5. 

5. 

Tsp 

TSP 

s 

Miscle  spindle  lead 
time  constant 

G(  1,7,8) 

0.02 

0.15 

1/11. 

1/11. 

TPS 

TPS 

s 

Hjscle  spindle  lag 
time  constant 

G(l,7,8) 

0.0 

0.1 

0.0 

0.0 

Tss 

TSS 

s 

Mjscle  spindle  hicfi 
frequency  lead  time 
constantly 

G(l,7,8) 

0.0 

0.1 

0.0 

0.0 

JKG 

N/N 

Golgi  tendon  organ 
model  gain 

G(l,7,8) 

' 

0.1 

1.0 

0.5 

— 

0.5 

T9 

IB 

s 

Golgi  tendon  organ 
time  constant 

G(l,7,8) 

0.02 

0.1 

1/18. 

1/18. 

h 

TZ 

s 

Golgi  tendon  organ 
time  constant 

G(l,7,8) 

0.0 

0.1 

0.0 

0.0  | 

[ 

TP 

TP 

s 

Golgi  tendon  organ 
lag  time  constant 

G(  1,7,8) 

0.0 

_ 

0.1 

0.0 

0.0  | 

1 

*Sources  Legend:  Measired  (e.g.,  photo)  or  specified  by  situation 
Tabulated,  cadaver  or  anthropometric  tables 
"Fitted,  to  transfer  firction  data 
"Sues st irate  (on  basis  of  physical  properties) 
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TABLE  1.  (Concluded) 


i  i  i - r 

— r 

- r 

— israwai — r 

- TYPICAL  VALUES 

~T 

iiii 

1 

i 

RAAGE  | 

(STAfOARD  CASES) 

l 

1  !  1  1 

I  PARA-  |  FORTH*  |  UNITS  | 

LOCATION 

1 

1 

r 

SOURCE*  | 

1  1 

1  1 

1 

1 

r 

i 

1  PETER  I  NAPE  |  | 

IIII 

Definition 

1 

1 

(Refs.)  | 

1 

MNIMJM  |  HVXIMJM  | 

1  1 

sum  i  sidcrw 
l 

i 

i 

1 

1 

1 

1 

1 

1 

1 

1  IMAGE  FIXATIONS  STI BULO- 
1 - OCUL*  SERVO - 

1 

1 

I 

1 

1 

1 

1 

r 

i 

i 

i 

1 

1 

1 

r 

i 

i 

i 

1  KOE 

! 

i 

i 

i 

| 

MODE 

1  rad 

1  ra? 

1 

I 

1 

I  Fixation  error  gain 

1 

1 

1 

1  F  (9) 

1 

1 

1 

1 

1 

1 

1 

1 

1.0  | 

1 

1 

1 

10.0 

1 

1 

1 

1 

7.06 

I  7.06 

i 

i 

i 

i 

i 

j  *RE 

1 

i 

i 

i 

i 

XKRE 

1 

1 

I  rad 

1  rails 

1 

1 

I  Fixation  rate  gain 

1 

i 

1 

1  F  (9) 

1 

I 

1 

1 

1 

1 

1 

1 

0.0  | 

1 

1 

10.0 

1 

1 

1 

1 

1 

0.0 

I  0.0 

1 

|  i 

i 

i 

i 

[  TV 

I 

i 

I 

i 

TV 

1 

1  s 

1 

1 

i 

(  Tims  da  lay  in  fixation 

1  loop 

1 

1  F  (9) 

1 

1 

1 

1 

1 

0.0  | 

1 

1 

0.10 

1 

1 

1 

1 

0.045 

I  0.045 

1 

i 

i  i 

1  a 

I 

i 

i 

ALPHA 

1  rad/s 

1 

| 

I  "Pseudo  Integrator"  bre* 

1  frequency 

1 

1  G,  F(9) 

1 

| 

1 

1 

l 

0.0  ( 

1 

1 

0.5 

1 

1 

1 

0.3 

1  0.3 

i 

i 

l  i 

|s 

l 

I 

i 

I 

XKP 

I  rad 

I  Fa? 

i 

1  Ta-get  pirsuit  gain 

1  (Optional) 

1 

1  G{9) 

1 

1 

1 

1 

1 

1 

0.0  | 

1 

1 

1.0 

1 

1 

1 

1 

0.0 

1  0.0 

1  j 

I  *IE 

r 

i 

i 

i 

HUE 

1 

1 

I 

I  Position  gain  frm 

1  "Vestibular"  sensor 
| 

1  F(9) 

1 

1 

~r~ 

i 

i 

1 

1 

0.1  | 

1 

1.0 

T 

1 

1 

1 

0.67 

1  0.67 

Jf 

j  i 

r* 

i 

i 

| 

MOVE 

1 

1 

1 

1 

1  Velocity  gain  from 

I  "  Vestibule"  sensor 

1 

|  F  (9) 

1 

1 

i 

1 

0.0  | 
j 

0.1 

1 

1 

1 

0.1 

j  °'1 

! 

i 

1  tc2 

i 

i 

TC2 

1 

1 

1 

i 

I  "Vestibula'"  lag 

1  F,  T(9) 

1 

1 

i 

_i _ 

0.0051 

1 

0.015 

1 

1 

1 

o.m 

I  0.01 

i 

i 

i 

r 

i 

“TW 

1  s 
| 

|  Ocular  Serw  Lag  time 
constant 

rn?) 

| 

| 

(ojrr 

1 

nr 

T 

1 

in - 

I - O — 

| 

i 

1 

J  tlm 

i 

i 

i 

TLM 

1  s 

1 

1 

1  Ocular  Servo  Lead  time 

1  constant 

1  F  (9) 

1 

1 

i 

i 

i 

0.0101 

1 

1 

0.020 

1 

1 

1 

1 

0.0125 

1  0.0125 

1 

1 

1 

1 

i 

i 

!<* 

1  ¥ 

1 

1 

1 

1 

j 

1 

1 

1 

Switch  for  fast  mode  in 
Ocular  Servo 

1.  use  S,l€ 

0.  VE  *  infinity 

1  M 

1 

1 

1 

1 

I  0.0 

1  (V00) 

1 

1 

1 

1 

1 

1 

1.0 

I  1.0 

1 

1 

1 

| 

1  1.0 

1 

1 

1 

| 

1  1 

1 

1 

1 

I 

1 

| 

jce 

1 

1  2E 

| 

Ocular  Servo  Danping  Ratio 

1  F  (9) 

1 

I  0.2 

1 

1 

1 

1.0 

|  0.65 
| 

I  0.65 

| 

1 

1 

j 

1  u 

1  o 

1 

1  V£ 

1 

|  rad/s 

|  Ocular  Servo  Natural 

1  F  (9) 

1  200.0 

1 

1 

400.0 

1316.23 

|  316.23 

1 

I 

■  ® 

1 

1 

1 

1 

frequency 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

J 

♦Soirees  Legend:  Iteasured  (e.g.,  photo)  or  specified  by  situ*1on 
T*ulated,  ca<hver  or  anthrepcnetrlc  tables 
Titted,  to  transfer  function  data 
‘Guesstimate  (on  basis  of  physical  properties) 
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TABLE  2. 

PARAMETER 

FILE 

STRUCTURE 

(SHORT  FORMAT) 

A. 

Lower  Body 

Mb 

Bb 

kb 

kbs 

eL 

G 

B. 

Torso 

Mp 

Btb 

Kps 

®T 

% 

0TN 

%N 

9s 

Bs 

C. 

Neck 

% 

In 

but 

Knt 

®N 

% 

D. 

Head/Display  Viewing 

% 

% 

®HN 

kHN 

cH 

% 

By 

VD 

E. 

Arms  (Upper  then  lower) 

h 

ei 

L1 

M2 

12 

9g 

L2 

W1^ 

F. 

Grip  Interface 

eU 

bi/ki 

CI 

G. 

Stick 

Ms 

BS 

KS 

CKS 

Sx 

sz 

ec 

Lc 

Ksc 
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APPENDIX  B 

DEFINITION  OF  FORCING  FUNCTION  AND 
RESPONSE  VARIABLES 


TABLE  3  INPUT  (FORCING  FUNCTION)  VARIABLES 


1 

|  MATRIX 
|  COLUMN 
|  NUMBER 

1 

1 

INPUT 

VARIABLE 

NAME 

COLUMN 

CODE* 

UNITS 

POSITIVE 

DIRECTION 

DEFINITION  AND  USAGE 

1 

1 

1 

1  1 

1 

Nzp 

NZP 

g 

Upward 

Platform  acceleration 

1  2 

1 

DZP 

m 

Upward 

Platform  deflection 

1 

1  3 

1 

Nxp 

NXP 

g 

Forward 

Platform  acceleration 

1 

1  4 

1 

% 

DXP 

m 

Forward 

Platform  deflection 

1 

1  5 

1 

0FD 

N 

Forward 

Force  input  to  stick 

1 

1  6 

1 

Bzp 

BZP 

m 

Upward 

Special  Platform  deflection* 

1 

1  7 

| 

Bxp 

BXP 

m 

Forward 

Special  Platform  deflection* 

1 

1  8 

NMC 

NMC 

N 

Neuromuscular  command  within 

1 

CNS 

1  9 

1 

1 

9ti 

TTI 

rad 

Test  input  into  fixation  error 

♦Small  "b"  represents  a  typed-blank  space. 

Inertial  shoulder  (or  head)  motion 

♦These  are  used  to  get  transfer  functions  o7  Platform  motion 
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TABLE  4  OUTPUT  (RESPONSE)  VARIABLES 


MATRIX 

OOliMI 

NLMBER 

RESPOND 

VARIABLE 

WE 

wvmix 

CCLUW 

CODE* 

UNITS 

POSITION 

DIRECTION 

DEFINITION  AND  USAGE 

1 

c 

f£|t 

m  or  N 

Forward 

Stick  output  (Kfs  used  to  scale  as  foroe) 

2 

^ca 

FCA 

N 

Forward  on 

Interface  foroe  at  stick  yip 

. 

stick 

3 

Aq 

pa 

m 

Upward 

Vertical  shoulds”  deflection  relative  to 

platform  (zs-Zp) 

4 

A01  + 

an. 

rad 

Pitch  ip 

Lows’'  ami  angle 

5 

A0a 

OTA 

rad 

Pitch  ip 

Lows'  arm  angle 

6 

fly 

F1Y 

rad 

Ip  on 

Vertical  interface  force  at  elbow 

upper  arm 

7 

N 

Muscle  nrodel  foroe  referred  to 

upper  arm  c.g. 

8 

Li 

*LI 

N 

Forward 

Grp  interface  deflection  (stick  to 

lower  am) 

9 

fix 

FIX 

N 

Forward  <n 

Horizontal  interface  foroe  at  elbow 

upper  arm 

10 

% 

|A)M 

m 

Internal  state  in  muscle  model 

11 

fn 

FNS 

N 

Forward  on 

Horizontal  interface  foroe  at  shoulder/ 

ipper  am 

ipper  am 

12 

fvs 

FVS 

N 

Ip  on  ipper 

Vertical  interface  foroe  at  shoulder/ 

arm 

ipper  am 

13 

*s 

0XS 

m 

Forward 

Horizontal  shouldr  deflection  relative 

to  platform  (xj-Xp) 

14 

a&t 

an 

rad 

Pitch  ip 

Torso  rotation 

15 

1 

ffh 

Ftti 

N 

Forward  cn 

Head/heck  interface  force  | 

)  16 

1 

A6fl 

OTH 

rad 

Pitch  up 

1 

Head  rotation  I 

1  17 
| 

fyh 

FVH 

N 

Up  cn  head 

Head/neck  interface  force  | 

J 

1  18 

1 

fby 

FBY 

N 

Up  cn  Torso 

1 

Vertical  hips/torso  interface  force  1 

1 

1  19 

't* 

FBX 

N 

i  i 

Forward  on  |  Hips/torso  Interface  force  1 

1 _ 

hips 

1 

>  *ttote:  Snail  V  represents  a  typed -bl art;  space. 
iNote:  Afy  (Colim  4)  naans  snail  signal  perturbation  cf  0\  given  In  fig.  3 
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TABLE  4  (Continued) 


MATRIX 

RESPONSE 

MATRIX 

POSITIVE 

COLUMN 

VARIABLE 

(DlifN 

UNITS 

DIRECTION 

DEFINITION  AND  USAGE 

NUMBER 

NW€ 

OOOE* 

20 

“5p 

UEP 

m 

Up  the  seat 
back 

Internal  state  in  hip  nodal 

21 

‘tp 

UBP 

m 

Up  the  seat 

Hip  deflection  relative  to 

back 

platform  (parallel  to  seat- 
back) 

22 

xi 

(®CL 

m 

Internal  state  In  nuscle 
nodel 

23 

a^N 

m 

rad 

Pito)  ip 

Neck  rotation 

24 

fnN 

FNN 

N 

Forward  on 

Neck/torso  horizontal  inter- 

neck 

face  force 

25 

fVN 

FVN 

N 

Up  on  neck 

Neck/torso  vertical  inter¬ 
face  force 

26 

Rrt) 

RW) 

tfeg 

Up  relative 

Head  point  of  Regard  relative 

to  displey 

to  a  display  on  the  platform  | 

Z 7 

(WA 

N 

Muscle  nodel  caimand  out  of  | 
spinal  cord  | 

28 

fa 

IdFA 

N 

Commanded  force  out  of  force  | 
activation  dynamics  | 

29 

Faa 

FAA 

N 

Internal  state  in  force  activa-j 
tion  | 

30 

*9 

J4AG 

N 

Output  of  gplgi  tendon  organ  | 

sensors  | 

31 

* 

VfiC 

N 

Spinal  cord  ccmrend  | 

32 

Ae 

N 

Error  between  neu-onuscula*  | 

caimand  and  proprioceptive  | 

stick  force  | 

33 

xsp 

XSP 

m 

Muscle  length  change  sensed  | 

by  spindle  nodel  I 

34 

Arn 

DRN 

m 

Relative  normal  deflection  of  | 
am  and  rest  I 

35 

Art 

CRT 

m 

Relative  tangential  deflec-  | 

tion  of  am  and  rest  | 

36 

A2h 

D3i 

m 

Vertical  conponent  of  head  | 

notion  relative  to  platform  | 

♦Note:  Small  *4"  r presents  a  typed -blank  space. 
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'C*  1T» 


/. \jZ  . 


TABLE  4  (Concluded) 


MATRIX 

COLUm 

NUMBER 

RESPONSE 

VARIABLE 

NAME 

MATRIX 

COLUMN 

CODE* 

UNITS 

POSITIVE 

DIRECTION 

DEFINITION  AND  USAGE 

37 

*cl 

AC1 

rad 

Output  of  compensatory  action  in 
eye  fixation  loop 

38 

*E1 

AEI 

rad 

Output  of  "Pseudo  Integrator" 

39 

^MC 

TMC 

rad 

Coirmand  to  Ocular  Servo 

40 

Ay 

bAV 

r/s 

Head  up 

Output  of  "Vestibular"  Sensor 

41 

eEH 

TEH 

rad 

Up 

Eye  rotation  relative  to  head 

42 

6M1 

TM1 

rad 

Internal  node  in  Ocular  Servo 

43 

e. 

THI 

rad 

Up 

Fixation  error 

44 

RED 

RED 

m 

Eye  up  rel. 
to  display 

Relative  eyem  point  of  regard  at 
display 

45 

eTC 

TTC 

rad 

Output  of  "Image  Pursuit"  opera¬ 
tion 

46 

ZT 

HZT 

m 

Up 

Display  inertial  displacement 
(vertical) 

47 

eE 

THE 

rad 

Up 

Inertial  eye  rotation 

48 

DXH 

m 

Forward 

Horizontal  component  of  head 
motion  relative  to  platform 

♦Note:  Small  "tf"  represents  a  typed-blank  space. 


APPENDIX  C 
EQUATIONS  OF  MOTION 

if  cKs  4  0 

(sin  ec)c  -  xs  -  L1(cose1)-^e1  +  L2(sinea)Aea  -  LjCcobGjj)  =  -  Bjjp 

if  cKs  40  it  CKs  =  0 

f - 7* - — — -  — 

|-(MgS2  +  Bss  +  Kg)  -  ^  -  M2(cosec)  s2  J  c/KgC  -  c 

+  (sinec)Lca  +  (cos0c)lr|y  -  M2(cosec) (sineIj)s2LI 
-  M2%(cos6c)  (cosea)s2^ea  +  (cosGcKB^s  +  KAR)^n) 


+  (cos0c)(Bfs  +  KF)(Art) 

=  “  (  2/s  +  M2  coS0c)aZp  -  Sx/g  -  Fd 

3-  • 

( CO 80  ) 

Az1  +  L|  (sine1  )ZiQ1  +  ^(cose^^  +  Ljisinejj)  +  — c  *  0 
k  |  I^2  +  (%M-,g  +  L-jM^  -  I^Favo)0080!  |  -  L^sinG-,)^ 


+  Lb(coss1  )f^  +  La(cos01  )fns  +  La(sin@1  )£vs  "  siLaFM  -  0 


NOTE:  Favo  =  Arn^  (M,  +  M2)g  \  0a  =  90  +  0!  -  6ej 

**p  =  SNZp  +  s2(^p  +  BSp)  >  =  gN^t  +  s2(AXp  +  B^) 
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5 


j(l2  +  M2l^cos2ea)s2  +  (FaV()  -  M2g)Lesin0a  -  I^fsinejf^ 


-  L2(cosoa)f1y  -  Ld( sinea)fCa 

~  +  Ldcos20a)(BARs  +  Ky^)^rn 

-  (l^sin0acos0a)(BFs  +  KF)£rT 

if  cKg  ^  0 

+  IflM2(cos0a) (cos0c)s2c  +  %(cos0a)M2(sin0jj)s2Lj  =  cosea)M2az 

p 

f1  +  fVs  ~  M1s2Az1  -  M1La(sine1)s2^01  =  Mja- 

Jr  * 


Cj(cosQjj)fc^  +  |cj( sine jj) %2s2  +  Bj/Kj  +  1  ^  Lj 
+  CI(sin0IJ)(cos0a)M2Lds2^a 

-  CjCsinejj)^ 

- iW - 

+  CjCsiiiCjj)  [-  (Bfs  +  KF)(sinea)Art] 

=  CjCsinOj^Mgazp 

+  CjfsinOjj)  [-  (B^s  +  KAR)(cosea)Arn] 


if  Ck_  ^  0 

4-  CI(sin0IJ)  (cos0c)M2s2c 


9  if  cKs  ^  0 

-M2(sinec)s2c  +  M2(coseIJ)  s2L-j-  -  |  (M?Ldsin0a)s2  +  FaVo  |  Aea 
^sc 

+  fix  -  fCa  “  (BFs  +  KF)(cos0a)£rt  +  (B^s  +  K^)  (sinea)£rn  =  M2^ 
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fns  “  f1x  ~  Ml  s^xs  "  M1La(cos81)s2^e1  = 

fvs  +  f VN  ~  fby  +  Mr(coseL)s2Ubp  -  MrLT(sineT)s2^6T  =  “  MraZp 
xs  +  (sineL)ubp  +  (Lpcosej  +  Lscoses)A0T  = 

[LpS2  +  (Bpg  +  Bjg'jOs  +  10pg  +  Kot  —  WE]  A0ij,  —  (BjuS  +  Kjjiji)Aejj 

—  (Lscos8s)fns  —  cos  9^)1*^  —  Ls(sin0s)fvs 

-  Lrflj(sin&jTpj )fyjj  +  Lj.(cos9r)fbx 

—  I»ji(sin&p)fby  +  S^IiaFj^  =  0 

WE  =  lTC°s0rf(%  +  J%  +  Mj,  +  M-|  +  MgJe  —  FaVQ] 

+  Lscos0s[(Mt  +  Mg)g  —  FSVq]  +  IrjjcosajTjj(^{  +  Mjy)g 

fnb  +  Mg^SinOg)^  J^jp  +  Mgr(LpCOSft]i  +  I/pjjCOS0r]ijj)s?A0kji 
+  cos0jj)s2A9jj  +  MjjLg(cos0g)s2‘Ci8g  = 

if  NM  0 

[(Ihs2  +  Bhjjs  -  MggLHCOSOH)^  +  KHN]^H_+_ * 

-(BhnChs  +  kHn)A0N  ~  Ch^C3100!!)^  ~  =  0 

fvh  —  Mfj(  0030^)32^^  +  %(Lipsin0j  +  L^sin9j^)s2^0j 

+  %LN(sin%)s2A0N  +  MfiLj|(sin0ji)s2-^0g  »  ^arp 
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l8 


19 


(coseL)fby  +  (Bbs  +  Kg^p  +  (sineL)fbx  +  MgS2^  =  Mg^ix^)^ 

-  MgCcoseL)^ 

+  fnN  +  fns  -  cos  0^32/^,  _  Mr^sineL^Ubp  =  -  Mj.aXp 


20 


/  BflS  +  Kb 

\  kBs 


U5P  “  ’  bP 


0 


21 


(coseL)ubp  -  (Lrp3in0r  +  Lssin0g)A0T  -  Az-| 


23 

[%s2  +  b1JTs  +  KNT  ~  (^^8%  +  )cos0n]a9n 
~  (ainen)fvN  ~  %l(cos0N)fnN  "  “  Iwi)sin0N  +  I^sinon] 

-  [LjjCosBjj  +  (%  -  LN^cosOjjjfjjh  +  [lHs2  - 

(BjjrpS  +  Kjqrji)^0rp  =  0 

24 

fnN  —  *nh  +  *%(sin0L)s2ubp  +  Mjj(I^cos0j  +  X/jijjcos0^jj)s2Aft|, 

+  MN(LHlcos0N)s2A0N  =  %axp 


25 


fyjj  —  ^vh  —  ^%(cos0L)s2ub  +  I%(Brsin9j  +  Ljjj3in9pjj)s2A0fp 
+  %(I*J1sin%)s2AeN  =  %RZp 
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NEUROMUSCULAR  MODEL  EQUATIONS 


7 


Fm  +  (BrpS  +  KT)La(Ae1)  -  (ftj.8  +  Kp)%  =  0 


10 


Fm  +  (Bes  +  K^Xm  -  (Bes  + 


0 


22 


Fm  -  (Bms  +  KM)La(A9T)  +  (Bms  +  %)xt  +  FA  =  0 

if  NM  =  0 


if  NM  =  0 


27 


(Tpss  +  1  )(Aa  +  Ag  —  Ac)  ~  KSp(TSpS  +  1  )(Tsss  +  1  )xSp 


28 


(xas  +  2)FA  +  Ka(Tas  -  2)Faa  =  0 


29 


F  —  F 
aa  raa 


This  set  to  1.0 
if  =  0 


if  NM  =  0 


50  (Tps  +  1)A  +  Kg(Tgs  +  l)(Tzs  +  1 )%  =  0 


51 


(tcss  +  2)(Tccs  +  1  )AC  +  K1c(tcss  2)Ag  0 
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-%<W  *1*0* 


32 


NM  =  0 

/_ cos'  nc s‘’  |  c  —  slnpcfca  +  (M2I^coseacosecs2)Aea  -  cosoc  fiy 
VKsc  ' 

o 

+  (MosinO|;jcosocs‘'')LI  -  (B^coseaCosepS  +  KARcoseacosec)'!Vn 
-  (BFsinoacos0cs  +  KpsineacosPc)Art  +  Ap 
=  M2CosecaZp  +  NMC 
NM  =  1 

Ae  +  RHD/VT)  =  NMC 


33 


xSp  +  ~  ^  =  0 


26 


[I^sinfoT  —  Py)  +  l*j^jsin(0'pj^  —  Qy^^Of  —  Vjy^GH  cos(0E  —  0y)u)jp 


+  LNSin(oN  -  Oy)^N  +  RHD  =  BZpcosOy  -  B^sinGy 


3** 


A r  + 

^  n 


if  cKs  f  0 
c  "cos(e  ~-^0a) 
^sc 


+  (L2  -  3JER)^Ga  +  Lisin(ea  -  ejj)  =  0 


3l> 


Art  -  c 


if  cKs  4  0 
^sc 


+  LjCOs( Ga  -  Bjj)  = 
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36 


+  ( Xjpsin0ip  ^ 


-  (cose^Ujjp  +  (%sineN)^eN  =  Bgp 

57  (2  +  V5)^-!  +  (Tvs  -  2)(%es  +  %E)eI  =  0 


38 

39 

1+0 

1+1 

1+2 

1+3 


(s  +  a)AEI  +  A  +  KjeAv  =  0 


eMC  ~  0TC  “  AEI  +  %EAV  =  0 


(tC2s  +  1  )fiv  -  s^eH  =  0 


Qe  =  0  |  0EH  "  eM,  =  0 

%  =  1  {  ^  +  S  +  ^6EH  "  0M1  =  ° 

e  ' 


(1  +  TEMs)8Ml  -  (1  +  ^8)0^  =  0 


RED 

01  ~  vj 


®TI ' 


1+1+ 


[LjSinte-r  -  0y)  +  L^sinCe^  -  0y)]^T  -  VpAeji 


-vkbpcos(eL  -  &y)  +  AeuLnsinteN  -  ©y)  -  vDeEH  + 


Bz  cosey  -  BxpSinev 


RED 
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U5 


Kp 

®TC  _  ZD 


0 
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eE  “  ~  °EH 


0 
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Axjj  +  (sineL)ubp  +  (Lj,coseT  +  I/j^cose^^T  +  l^(coseN)^eN 
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APPENDIX  D 


TYMSHARE  EXAMPLE  PROBLEM 


This  appendix  documents  the  use  of  CREATE,  BIODYN  and  PLOT  on  the 
Tymshare  System  31  PDP  10.*  The  computer  dialog  appears  on  the  next 
several  pages  for  a  typical  session.  An  existing  PARAMETER  file  is 
modified,  a  CHOICES  file  is  created,  BIODYN  generates  the  requested 
transfer  functions  (only  one  of  which  is  illustrated),  and  PLOT  produces 
a  quick-look  Bode  plot.  Throughout  the  dialog,  all  user  inputs/  respon¬ 
ses  are  underlined.  The  following  6  character  filenames  are  accessed  in 
the  course  of  this  example  (these  are  not  the  Standard  Pilot  and  Stan¬ 
dard  Crewman  which  are  available  on  Tymshare): 


File  Name  Title 


Existing  PARAMETER  file: 
Modified  PARAMETER  file: 
CHOICES  file: 

TF  file: 


BB19 

"Stiff  Stick 

SEMI 

"Semi-Supine 

CHOSAR 

TAPE19 

There  are  several  differences  between  the  CDC  version  of  the  BIODYN- 
80  package  and  its  Tymshare  counterpart.  These  are  listed  below: 

•  All  three  programs  run  Interactively  on  Tymshare;  all 
files  are  specified  by  the  user  during  execution. 

•  CHOICES  filename  can  be  an  existing  file,  whose  name 
is  to  be  modified.  (In  CDC,  the  CHOICES  filename 
must  be  a  new  name.) 


*A  user  planning  to  use  BIODYN-80  on  Tymshare  should  be  familiar 
with  the  Tymshare  manual,  XEXEC,  especially  Section  3.  At  present, 
BIODYN-80  is  not  a  current  Tymshare  Library  Program,  and  a  potential 
user  ahould  contact  the  second  author  for  the  required  procedures. 


•  Because  of  the  above  differences,  there  are  slight 
changes  In  the  dialog  for  the  Tymshare  version  (e.g., 
the  NEW  PARAMETER  PILE?  query  is  not  used,  since  the 
PDP-10  software  can  determine  internally  whether  or 
not  a  new  filename  has  been  input). 
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D  •  1  Running  CREATE 


+ 1 


-LXL  LIGATE 


LEAPING 
'"Yf  CUT  ION 


TNri.iT  PARAMETER  FILENAME 


tl<  19 


LISTING  DESIRED  ? 

(LO>NG  OR  t SH ) OR  T  LISTING' 


M 

PARAMETERS  F OR  I 
STIFF  STICK 
LOWER  BODY 


11.000  2290.1 

TORSO 

29595. 

?  1  r / 1  ? . 

13.000 

18.000  0.80000 

! 6 . rso 

500 . 00 

10.000  0  •  1 11000 

5.0000 

0 . 30000 

5 . 0000 

NECK 

0.00000  0.00000 

0.25000 

00 . 000 

■20.000  0.10000 
HEAP/ DISPLAY  VIEWING 

0.50000C  01 

1.3400  0.39000F *01 

0.00000 

l  .  0000 

o.ooooo 

70.000  0.00000 

ARM  (UR'  ER,  LUWLK  > 

-30 . 000 

0 . 70000 

1.3720  0^120001-01 

10.000 

0 . 290 AO 

0  11000 

1.017'  0*.  157O0E  01 

i  v*.  .oo 

0 . 30500 

0.50000 

GRIT  INTERFACE 

0.0000  -  0. 1  C'OOOfc  -01 

0. 5579  IE -01 

•jT  icf 

0.11000  2.0)00 

13900. 

l .0000 

0.00000 

90.00“  0. A 1000 

ARM  F\ !“  "i  T 

13900. 

0.00000  O.'MOOO 

0 . 00000 

0.00000 

0,00000 

/ it  UkGhi , :•  r Ul  AP  .  f  S  T FM 

A  01. 40 

0.2013  V -01  •».90709E*01 

0.R938RE  -01 

1.0000  •••.00000 

14.000 

0 . 00000 

1  ,0A?0  HOO 

7.4312 

40.000 

0.00000 

1*1. 0000  0.'.  OVOOf-OI 

A . 0O000 

0 . 00000 

0,50)00  H.5G554F-01 

0,00000 

0.00000 

IMAGE  i  rXATTOrJ/vesriPHi 

i.) -O'  HI  AH'  SERVO 

7.0.*vr.  O.'OOOO 

0.  15000 E-01 

>.  30000 

•'.00000 

«>.  A  7000  0.10000 

0. 1O000E  o 1 

o.ioo  :  .oot -<  i 

1  .uOOO 

•  • . 65000 

<14.23 

FHAflGFS  PL'  -  IRFTl  - 

yr 

[Mr  in  MNFMGNIf  *j  CHARS)  AND  ‘-'All 

!|"  FOR  f  ACM 

CHANGE f 

r«)  '  x  t  T »  t  yf  r  xxx. 

9.8000 

0.30000 


0.00000 

0.00000 


Execute  CREATE 


Interactive  input  of 
PARAMETER  file  to  be 
modified 


y  Generate  short  listing 

■f  of  file  to  modify 


Make  changes  to  file 


no. 


(continued  on  next  page) 
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D.1  (Continued) 


■  Ui.Vi.  i.rsiKM-' 


a  ,i'M,  tm  I’jH-iihr  i.  ;srirn_',? 

H 

■  .ms  ii  i  ■ ,  r  cvc 

:  ill  -  '!if-  i Nt. 


j-i  r pm 


:  .  ■••■o 

IT'-'O.  i 

'LC19. 

1 

u  .  «>co 

0.0(000 

\  • .  *’V. 

•00.00 

aZ .050 

r:  ■ 

. .0000 

0 .  50000 

. , ‘ '00 

■:» .  ooo  *;• 

.1.0  0  7 

■  0../00 

'Vj 

0.  ICO''' 

C  .  ■  00*.  ■!:  *01 

E  nl  M*’M  .V 

0  1 1 K  i  10 

; .  : ; ;  0 

0  •  .*•  9  00  Oh  -0  1 

■> .  OvOOO 

1 .0000 

o .  :.oooo 

i :  •■> .  oo 

o .  O’**  >oo 

:  >.000 

0. 75000 

Mi  i.TFM'. 

I  '  i\ 

•  . 

.  i 'OCT  •  «)1 

’•.000 

0 . J9000 

0. 1  V 

1  .  ;  1  "  ’• 

.  i  ..  m 

:  r.i.oo 

O.SO'/OO 

• 

:• :  :i .  r.  i 

r- 

i  •: 

o.  >  .  ’fur  oi 

.  «'■<>■*> 

.  C:  0  **  ? 

1  .*^00. 

1  .0000 

0 ♦ 00000 

0  .  M  ’ r 

:  W'-'O. 

r-i  -;t 

0.  . 

•  <>''-000 

0.00000 

0 . 00000 

■:  :r  )  '  ;ui 

!•  ’  .  :  1!  M 

i  1.  10 

.  '-UT  1 

0.  •><<■- ‘Wh  -01 

O.HVHBrtL  01 

}  .  jOO 

C 

1ft.  000 

v) .  00000 

o:;oc 

’.o 

• .  hi: 

10 . ooo 

0 . 00000 

"• .  0000 

<■ .  Vo  -<>* 

'  .00000 

0.00000 

.  ,  >'*00 

/  si  -  a 

.00000 

0.00000 

<*'  ! 

i,j/or  • :  ii.:l 

0  I.UJI.Y.R  SEIVOO 

'  .<■)<><,  oo 

■/.  no (>•:;•  f  oi. 

0. JOOOO 

0 . 00 ooo 

.  '•  <>0t 

0  .  ( ■»<  -  > 

'. 1 0000 H-Ol 

*/.  1  •••■  ,,%l  01 

l .0000 

•).  A7.000 

51ft. JJ 

r :  r  .*»mi 

f-jk  r.,HOhK 

ri  a  f  1 1  ie.  :  » 

Modify  title  for 
PARAMETER  file 


Short  listing  of 
modified  file 


Interactive  input 
of  filename  on  which 
to  store  modified 
PARAMETER  file 


(Continued  on  next  page) 
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D.1  (Concluded) 


INPUT  CHOICES  FILENAME 
CHQSAR 


BIODYN  TFS  DESIRED  FDR  FIVIB 


jjsl 


'Y 


transfer  function  input  : 

FIRST  LINE-RESPONSE  MNEMONIC  *  FORCING  FUNCTION  MNEMONIC 
<  AAA  * AAA )  i  ENTER  XXX  TO  STOP 
SECOND  LINE-PLOTTING  INFORMATION*  5  ITEMS  : 

BODE  LOWER  FREQ.  LIMIT 
BODE  UPPER  FREQ  *  LIMIT 
BODE  UPPER  PHASE  LIMIT  <0. 

BODE  LOWER  PHASE  LIMIT  <0. 

LIST  <1.  TO  LIST  TABLE*  0. 

IF  m  PLOT  DESIRED*  fcNTER  0. 

1 m- 


DEFAULTS  TO  200.) 
DEFAULTS  TO  -400.) 
FOR  NO  LIST) 

FOR  ALL  ITEMS 


1 . *  100. ,0. 

1.  f 100. *0. *0. *0. 

filtUl'If. 

i..ioo.»o.»o.»i. 


MA. 


Interactive  input  of 
CHOICES  filename 

No  PIVTB  TF's  desired 


Format  for  CHOICES  file 
information  input 


Input  CHOICES  file 
information 


EXIT 


Exit  CREATE 


D.2  Running  BIODYN 


1 


EXE  PDN80>N7RFN 

r*mm - 

ExE  CUT  l  ON 


Execute  BIODYN 


ENTER  INPUT  FILE  NAME :  SEMI 

enter  choices  file  name:  chqsar 

1 1  -M.ar-80  ll!53 


case:  semi -sup i me 

ENTER  SYSTEM  OUTPUTS  FILE  NAME'**  TATTOO 

nfu  tile 

F-  DAT  I  Nr.  UNDCRFIOU  FH  1II?4A 

FLOATING  UNDERFLOW  PC- 111? 46 


Interactive  input  of 
PARAMETER  and  CHOICES 
files 

BIODYN  execution 
Case  title 

Interactive  input  of 
TF  file 


denominator: 


« 1 C444E-26 


( 

16.591  > 

<  46.384 

)  (  56.079 

)  < 

97.512 

( 

100.00  ) 

(  285.86 

)  (  333.22 

) 

( ( 

. 44231 

8.3929 

>  3.7123 

7.5273 

)  ) 

<  < 

.35697 

1 1 .  .509 

»  2.9060  t 

10.929 

)  > 

<  ( 

.20726 

13.642 

t  2.8274  » 

13.346 

)  ) 

( < 

. 19065 

24.860 

t  4.9386  , 

24.365 

>  > 

( l 

.32015 

32.898 

*  10 . 532  t 

31 . 167 

)  ) 

<  ( 

» 83789 

32.931 

t  28.251 

16.921 

)  > 

<  < 

.19448 

57.685 

,  11.219 

56.583 

)  ) 

<  < 

. 65024 
. 18027E+11 

316.83 

»  206.02  r 

240.70 

)  ) 

numerator:  dth/dzp 


, 29296E* 27 

«  .28049 

<  . 30000  ) 

<  .00000 

>  (10. 

621  > 

16.590 

(  44.384  > 

<  97.067 

)  <100 

.00  ) 

-138.57 

<  209.05  ) 

<<  .44231 

»  8.3929 

3.7123 

*  7.5273 

)  ) 

(<  .21687 

»  11.488 

2 . 49  J  3 

,  11.214 

)  ) 

((  .98904 

.  16.317 

IA>»44 

,  7.39  16 

)  / 

'  <  .  us  t  30 

»  33.163 

28.232 

,  1/.400 

)  ) 

( (  . 20646  ♦ 

*  35.4  10 

11 . 440 

»  54.216 

)  ) 

( (  .65024 

»  316.83 

206.0? 

»  340.70 

)  ) 

\  (  •  »'4  4  5 

*  336 . 66 

.325 . 4 1 

.  86.371 

> ) 

.47996FF08 

.37719E 

•02 

mu  -123.31 

11-rt jr-80 

1  2  .*  58 

> 

> 


s 


EXIT 


Printout  of  DTH/DZP 
transfer  function 


(See  Figure  9  for 
explanation  of 
format ) 


Exit  BIODYN 
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D.3  Running  PLOT 


fc  \C  cnF'l.QTyCwUttt  CE*OPTF' 

0AM  NG 

F*n:iinON 

ENTER  TRFN  SYSTEM  TILE  NAME J  TAFE20 

rl  •  I  htF*  11-Mar- 80  11153 

&r. mi -  Slip  irit 

I'll 

c:  /NZf 
I l/DZP 
PTH/DZF’ 


Execute  PLOT 

Interactive  entry 
of  TF  file 
Case  title 

Dump  of  TF  file 
contents 


TAPE  jo 

,i ,  rhF-  i  i  Me  r-  no  1 1  :c;3 

i'  MI  ■SUP  INK 

UTM/DZP 
J»i  b 

r  1 1-:S T  'IM'ER  M’  01  1-5  CANCEL  I.Lt«  J 
16 . JV >000 
1 00*00000 


IRCONI’-'IKDER  I'tPOLl.S  CANlJfc L.L ED  i 
.  1 17  V3  VO0  *  0.5927*00 

.‘,502-%%00  *  316.83100 

Nl-MEPO !  01  f : 


PLOT  execution 


Transfer  function  DTH/DZP 

Identical  roots  in  numera¬ 
tor  and  denominator  are 
cancelled;  close  pairs  are 
retained 


.  JfXMV 
.05o0\>  ) 

4  . 00000 

>  (  -  10. 

671 

)  ( 

16.590 

v’.Od7  > 

< -138.57 

)  <  '09 

.05 

) 

>  ) 

(  . .?  1  A3  7 

[ 1 . 483 

,  2.491.4 

* 

11.214 

<(  , 9UVR4  » 

L/..3L7 

,  16.646 

» 

2.3917 

) ) 

/(  .05  MO 

J.  163 

,  28.237 

» 

17. 100 

>  ) 

(  .  .’*>64  A  • 

55.410 

,  11.110 

* 

54.216 

)  ) 

<•  , 'MVj  7  * 

136.46 

,  375.11 

* 

26.321 

1  ) 

■it  nom  i  nc.  i  nr> : 

•:  tA.vyt 

.  56.079 

>  <  v  ? 

513 

)  < 

285 . 86 

<  55  J.;  ’  > 

.  (  ,  _  ;  i  *,  ,J  7  • 

1 1 . JOV 

f  2.9060 

10.929 

>  ) 

j  (  . ; O' :  a  * 

l 3. 6 43 

t  2.8274 

13.346 

>  > 

(  \  ,  !  V.-'6’..>  * 

71,060 

,  4.9386 

* 

74.365 

1  ) 

<(  . ;:mj  * 

*2.098 

*  10.537 

* 

31.167 

)  > 

( •:  . Hj/nv 

5.’  .7  3  1 

»  79.251 

1 6 • 931 

)  ) 

<  (  .1/1411 

,  :  vi vt  o.' 

57.485 

f  1  1  ,217 

56.503 

>  > 

) 


> 


Transfer  function  to 
be  plotted  ( see 
Figure  1 1  for  explana 
tion  of  format) 


J 


(Continued  on  next  page) 
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D.3  (Continued) 


H.1MP  U-Hrfr-80  UJ53 

SEMI  SUPINE 


DTH/D2P 

V*t 


XXX  =  MAGNITUDE  (DEO 


40.  30. 


0.  -20.  -40,  -60.  -80. 


28 ,  >  X 

X 


X  0 

X  0 

X  0 

X  0 

X  0 
0  X 

0  X 

0  X 


7.1  & 


5.0  0 

0 

4.0  0 

0 

0.5 -  0 

2.8  0 

0 
0 

2.0  O 

n 

u 

1.4  0 

0.2 -  0 

O 

1.0  0 


200.  100. 


0.  -100.  200.  -300.  100. 


000  PHASE 


(Continued  on  next  page) 


Bode  plot  of  DTH/DZP 
XXX  is  the  magnitude 


Frequencies  which  are 
denoted  as  hertz  (HZ) 
or  radians /second 
(R/s)  are  exact  at 


respectively 


000  is  the  phase 
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D.3  (Concluded) 


T  AF'K  20 

u-Mor-flo  11:53 

:-T  hi -  SUPINE 


PTM/DZP 

bit 


h/S 

1.00 
I  .  12 
3  .  26 

1 .41 
l  .53 
1.  "8 
2.00 
2.24 
2 . 5 1 
2.82 
3.1* 
3.55 
3 .  98 
4.47 
5.01 

5.42 
4.31 
7.08 
f  .  94 
3.91 

10.00 
11.22 
12.59 
14.13 
i  5 . 85 
l  7.  73 
19.  V5 
.•7.39 

70 .  I  3 
1  i  .  47 

55.48 
*.V  .31 
4  4.*7 
■ .  r. ,  •  2 
',4.7  5 
v'.  10 
••:>.  79 
>7,43 
09.  13 
« 00 . 00 
xi  r 


r»* 

-48.34 

46.31 
44.27 

-42.21 

-40.15 

-38.07 

-35.96 

-23.83 

-31.67 

-29.46 

-27.20 

-74.88 

-22.48 

-19.97 

-17.34 

-14.56 

-11.58 

-8.38 

-4.91 

-1.17 

2.85 

7.30 

12.34 
16.45 
18.74*' 
70  /  72 

23.31 

26.34 

23.34 
77.92 

26.04 
23.40 
70.78 
17.17 
:  4.43 
I  2 . 60 
!  ->.82 
8.59 
4.37 
1 . 35 
7.51 


PEG 

175.9 

175.4 

174.9 

174.3 

173.5 

172.8 

171 .9 

170.8 

169.7 

168.4 
167.0 
1 65  *  3 

163.4 
161.2 

158.6 

1 55 . 5 

151.7 

147.1 

141.1 

133.2 

122.9 

ito.o 

89.2 

56.6 

24.4 

-J  .7 
-27.3 
-60  •  1 
-104.6 
-150.8 
-189.7 
-221.9 
-246.8 
-264.3 

775.8 
-  205 . 8 

799.4 
-31  1  .4 
-320.3 
-377.3 
-333.2 


Listing  of  frequency  (rad/sec), 
amplitude  (dB)  and  phase  (deg) 
at  20  evenly-spaced  increments 
per  decade  for  DTH/PZP 

(This  listing  is  only  generated 
if  XL  =  1 .0  in  the  CHOICES  file 
for  this  transfer  function) 


Exit  PLOT 
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APPENDIX  E 


INTERCOM  EXAMPLE  PROBLEM 


This  appendix  presents  a  typical  terminal  session  in  the  Intercom 
4.7  operating  system,  accessing  the  CSA  mainframe  at  Wright-Patterson 
AFB's  ASD  Computer  Center.  The  dialog  is  annotated  so  that  a  typical 
user  will  easily  understand  the  basic  sequences  of  parameter  entry  and 
Job  steps.  All  user  responses  are  underlined;  each  is  terminated  by  a 
carriage  return. 

This  particular  session  was  an  exercise  using  a  9  cm  viewing  dis¬ 
tance  and  the  standard  crewman  (STDCRW)  parameter  set.  It  was  one  of  a 
series  of  runs  which  attempted  to  determine  optimal  display  distance  for 
minimizing  image  motion  of  a  vertically  vibrating  crewman.  The  example 
is  carried  far  enough  for  the  potential  user  to  see  how  the  programs 
Interface  with  INTERCOM  job  control  and  file  management  commands. 

The  steps  followed  in  the  investigation  are  listed  below: 

1)  Log  in  to  Intercom. 

2)  Attach  the  PARAMETER  file  to  be  modified  and  name 
it  TAPE20. 

3)  Attach  the  CREATE  program,  called  EXECRT . 

4)  Run  EXECRT,  make  changes  to  existing  file,  assem¬ 
ble  new  CHOICES  file. 

5)  Assemble  batch  job  to  run  BIODYN  and  PLOT,  using 
the  two  PARAMETER  file  (old  and  modified)  and  the 
new  CHOICES  file. 

6)  Submit  batch  job  to  input  queue. 

7)  When  job  completed,  list  the  output  file. 

The  user  is  advised  to  retain  in  his  permanent  files  only  those  which  he 
wishes  to  use  in  the  future,  and  to  frequently  purge  his  directory  of 
unneeded  parameter  and  output  files  in  order  to  forestall  job  failures 
resulting  from  file  space  overload. 
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Typical  INTERCOM  Terminal  Session 
user  responses  underlined 


ASP  COMPUTER  CENTER  INTERCOM  5.0 
SYSTEM  CSA 

'  PATE  05/20/00  T  T ME  1 4 . 4/, .  3V  . 

PLEASE  LOGIN 

ENTER  3-DIGII  TERMINAL  T  Ti  — 

05/20/80  LOGGEP  IN  AI  14.44.52. 
WITH  USER  IP  CK 
EQUIP1 /POR  T  12/001 

COMMANP-  REQUEST. TAPE7.*FF 

COMMAND-  REQUEST. TAPES. *PF 


1 .  Log  in  to  INTERCOM 


COMMAND-  ATTACH .  T  APE20 .  STOCKUJ 
PF  CYCLE  NO.  -  001 
COMMAND-  ATTACH » EXECRT 


PFN  IS 
EXECRT 

PF  CYCLE  NO.  -  001 
COMMANP-  EXECRT 


b 


2.  Attach  files  to  I/O  units 


3.  Attach  EXECRT( CREATE)  program  file 


NO 

YE_ 

LI 

LO 

YE 

NO 


L.  Execute  EXECRT  program  module 
new  file  7  j _  5*  Select  existing  data  file  (STDCHW) 

LISTING  DESIRED?^,  List  data  file 

(LO)NG  SR  (SH)ORT  L  ISTTNG'i’  J _ _ 

< L 0 ) NG  OR  <SH)ORT  '  T-1: r  l  Nl,  I 


7*  Program  ejects  unrecognizable 


NEW  CHOICES  FILE 


b 


BIOOYN  rrs  DESIRED  FOR  PIVI 


"?h 


8.  Long  listing  option  selected 

9.  Affirmative  on  CHOICES  option 

10.  Negative  on  PIVIB  option 


TRANSFER  FUNCTION  INPUT  ! 

FIRST  LINE-RESPONSE  MNEMONIC,  FORCING  FUNCTION  MNEMONIC 
(GAA.AAA)  i  ENTER  XXX  TO  STOP 
SECOND  LINE-PI. OTTTNG  INEORMAIIUN,  5  ITEMS  ! 

BODE  LOWER  FREQ.  LIMIT 

BODE  UPPER  FREQ.  LIMIT 

BODE  UPPER  PHASE  LIMIT  <0.  DEFAULTS  TO  200.) 

BODE  LOWER  PHASE  LIMIT  (0.  DEFAULTS  10  -400.) 

LIST  (1.  TO  LIST  TABLE.  0.  FOR  NO  LIST) 

IF  NO  FT  OT  DESIRED.  ENTER  0.  FOP  All.  ITEMS 

~FTNP  not  permissible,  pi  ease  RPiNPurh11*  Illegal  parameter  selected 
RHP.PZP  1 

.  1  f  poo  > .  0 . .  0 . .  1 « 1" 


MAX  FREQUENCY  RANGE  IS  .1  DFCADE* 
PLEASE  r.T  INf'Ul  f  NT  IRE  FINE 


12.  Program  screens  for  out-of-limit 
values 


1 1 1  99 1  r.9 :  ’ 0 1  ’ 1  - - -  13.  Finally,  valid  input 

■jLU - ’it.  Exit  from  CREATE 

STOP 


CP  SI  COWlTS  EXECUTION 


Session  continues  with  batch  run 
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Sample  INTERCOM  Batch  run 


COMMAND-  EDITOR 


f—  1 


.  . CREATE. S 


ENTER  LINES 
S  AR .  CM150000  *  d  TCSA  . 


L8Q00 15. RIEDEL. (213)679-2281 


ATTACH. TAPE? . CHOICE 1 . 

ATTACH. TAPES. UIFU9, 

ArfACH.EXEBIO. 

ATTACH. EXEPLT. 

EXEB 10 ( TAPE A . OU  TPUT . TAPE  7 . TAPES  > TAPE 1 9 . TAPE21 ) . 


REWIND.  fZFElT'. 

EXEPLTC OUTPUT. TAPE19) . 
«EOR 


2 


.  .SAVE.GQFILF .NOSEQ 


•  .STORE .GOP ILE .LQ00015 

CT  ID*  L300015  PFN-GOFILE 
CT  CY -  001  00000123  WORDS.: 

*  .end 


COMMAND-  BATCH. uOFILE. INPUT. HERE 


3 


h 


COMMAND-  FILES 


—LOCAL  FILES  — 

•1EXECRT  *00:-'ILE  *INFUT 
HTAPE3  HTAPE20 
— REMOTE  EXECUTING  JOBS— 
SARCK02 
COMMAND- 

TASK  , TEL  *  TERMINATED 
TASK  EXIT  WITH  OUTSTANDING  10 
PC  =10664<4 
PS* 170000 
R0-003770 


$6UTFUT  *TAFE7 


Note:  Output  files  generated  by  BIQDYN,  vnen  listed, 
are  identical  in  format  to  printouts  shown  in 
Appendix  D. 


.  Invoke  INTERCOM  editor 


.  Create  batch 
job  stream 


.  Preserve  file  containing 
Job  stream 

.  Batch  job  stream  to 
input  queue 


Job  terminated  by  system 
crash 
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